Localized Matter and Geometry of the Dirac Field. 
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Within the framework of classical field theory, the connection between the Dirac field as the field 
of matter and the spacetime metric is discussed. Polarization structure of the Dirac field is shown to 
be rich enough to determine the spacetime metric locally and to explain the emergence of observed 
matter as localized waveforms. The localization of the waveforms is explained as the result of the 
local time slowdown and the Lorentz contraction as a dynamic re-shaping of the waveforms in the 
course of their acceleration. A definition of mass as a limiting curvature of the spinor-induced 
metric is proposed. A view of the vacuum as a uniformly distributed unit invariant density of the 
Dirac field with an explicitly preserved invariance of the light cone is brought forward. Qualitative 
explanation of the observed charge asymmetry as the consequence of the dynamics of localization 
is given. The classical pion field is obtained as a manifestation of stresses, mass and charge flux in 
localized waveforms of the Dirac field. Some implications of the finite size of colliding objects for 
high-energy processes are discussed. A hypothesis that known internal degrees of freedom are the 
local spacetime (angular) coordinates that have no precise counterparts in Riemannian geometry is 
proposed. 

PACS numbers: 11.27. +d; 11.30.Na 
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INTRODUCTION 



The purpose of this paper is to show that the clas- 
sical theory of the Dirac field, considered as a primary 
form of matter can explain those important properties of 
the observed matter which so far remain a mystery when 
viewed from the perspective of quantum field theory. In 
the first place, these properties are localization of ele- 
mentary objects and the origin of their mass and finite 
size. Another, not less intriguing question is the origin 
of the observed charge asymmetry of normal matter - we 
find only small, heavy, positively charged protons (nuclei) 
and light, negatively charged, poorly localized electrons 
as the only stable particles around us. It seems that a 
possibility to answer these big questions has been over- 
looked at the early stage of field theory. 

An idea that the fields i/j(x) of matter themselves can 
immediately define the metric tensor g^ix) was brought 
forward by Wigner 1| and Sakharov Q. From the 
physics perspective, this idea is extremely sound; coordi- 
nates can be measured only through positions and shapes 
of material bodies. [In quoted works, tensor indices of g^ v 
were due to the derivatives d fJi ip(x).] It appears that the 
Dirac field builds up the metric of spacetime without re- 
sorting to ad hoc derivatives and it does this in such a 
way that the time slows down in the domains of magnified 
invariant density. This observation alone leads to a natu- 
ral and startlingly elegant answer to these big questions. 
Merely in the spirit of Huygens principle, this fact leads 
to self-localization of the Dirac wave forms into small, 
heavy, positively charged objects of finite size while leav- 
ing a negatively charged fraction of matter in the form 
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of an agile substance surrounding these small and heavy 
objects. It also changes the image of the Dirac sea as the 
vacuum - a uniformly distributed unit invariant density 
TZ = 1 is identified with goo = 1/1Z 2 = 1 and replaces a 
continuum of oscillators with an unbound energy spec- 
trum. 

From the perspective of the present work, the Dirac 
field is important, not as a special representation of the 
Lorentz group, but as a field that accurately describes the 
hydrogen atom. Lagrangian formalism is not used and 
no symmetry is assumed a priori. The main focus is on 
the possibility of deriving the most important properties 
of observed stable matter and its motion starting from 
the basic properties of the Dirac field and its equation 
of motion. No significant attempt to develop a formal 
perturbation theory that could have dealt with the finite 
size of particles as the in- and out- states of the quantum 
scattering process has been made so far. Once localized 
wave forms are found they immediately can be used as a 
basis for second quantization and their fields can serve as 
Heisenberg operators. The best prospect of this study is 
connected with the possibility to bridge the gap between 
point-like particles of classical electrodynamics and the 
plane waves of the quantum theory of scattering. 

The logic of the present work can be outlined as fol- 
lows: 

The stage is set in Sec lII A[ beginning from a review 
of well-known properties of the bilinear forms of the 
Dirac field with emphasis on their purely algebraic ori- 
gin. These forms are empirically verified to be affine 
Lorentz tensors at a generic point and they are further 
used to build a quadruple of orthogonal Lorentz unit vec- 
tors (tetrad). The possibility of treating these unit vec- 
tors as the tangent vectors of the coordinate lines of a 
usual holonomic coordinate system and thus to define the 
Riemannian metric as a descendant of the Dirac field is 
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considered in Sec III Bl It appears that certain conditions 
of integrability should be met and that these conditions 
are controlled by the Dirac equation. A key observation 
that the matter-induced metric must be equivalent to a 
long range interaction is made. 

The rules of differential calculus for the Dirac field in 
curved spacetime are reviewed in Appendix mostly 
following V. Fock A greater generality than in Q 
is intentionally admitted - and there was no possibil- 
ity to truncate it later. Sec lIIII thoroughly investigates 
if various differential identities, derived from the Dirac 
equation, can be put in the form of tensor equations, 
thus being independent of a particular choice of coor- 
dinate system. For identities that involve the energy- 
momentum T b a , the conclusion is negative with the fol- 
lowing facts firmly established: (i) The normal covariant 
form of the energy-momentum conservation cannot be 
assembled when the coordinate system is normal, (ii) 
The tetrad components T h a of energy-momentum are not 
the invariants of tensors, (iii) Being formally translated 
into coordinate form, the identity of energy-momentum 
balance keeps an explicit dependence on tetrad vectors. 
It does not reproduce the equation for a geodesic line in 
a given metric background, (iv) An explicit expression 
for the force of gravity (inertia) is derived from the con- 
straint, which accounts for the interplay between scalar 
and pseudoscalar quantities in the course of a physical 
acceleration of a Dirac object. Only in a crude approx- 
imation of a point-like object is the standard metric ex- 
pression recovered. To account for the flux of momenta 
in spacelike directions (pressure) in a localized waveform, 
a new stress tensor, Pg", is introduced and studied in the 
same detail in Appendix|B] A connection with the theory 
of Nambu and Jona-Lasinio is traced. The wave equation 
for the pseudoscalar density (pion field) with the source 
that has the structure of the axial anomaly is derived in 
Section [B2]. 

An unremovable dependence on tetrad vectors 
prompted a detailed investigation (in Sec lIV[) of the ge- 
ometric properties of vector and axial currents and con- 
straints that affect the integrability of differential equa- 
tions for their lines. It is shown that for the on-mass-shell 
Dirac field, the timelike congruence of lines of the vector 
current always is normal so that there always exists a 
system of hypersurfaces of a constant time. The key re- 
sult reads as dt — IZdso , where 1Z = \fp is the invariant 
density of the Dirac field; dt and dso are the intervals 
of the world and proper time, respectively. It immedi- 
ately predicts a general trend of self-localization for the 
Dirac field and a Lorentz contraction of accelerated ele- 
mentary objects as a physical process. Investigation of 
constraints connected with identities for the axial current 
(which, having a source, determines a radial direction) 
brought about another result - the maximal curvature of 
the 2-d surface of constant radius cannot exceed m, the 
mass parameter in the Dirac equation. A set of relations 
that connect bending of coordinate lines with the distri- 
bution of the axial current is obtained. These relations 



prompt a strong parallel between the local dynamics of 
the Dirac field and systems of inertial navigation - lin- 
ear acceleration inevitably causes a precession and vice 
versa. Sec lIV Cl deals with the intuitively appealing (and 
possibly not realistic) case of normal radial coordinate, in 
which behavior of the angular coordinates can be studied 
analytically. 

With the metric explicitly depending on the field of 
matter the Dirac equation becomes nonlinear in a unique 
way, which leads to self-localization as an intrinsic prop- 
erty of the Dirac field. Different forms of this equation 
along with an analysis of individual terms are the subject 
of Sec|V] and Appendix [C] 

In Sec lVH a major conjecture regarding the nature of 
electric charge is made. Maxwell equations are intro- 
duced and it is shown that a stable Dirac waveform can- 
not interact with its own electric field. Furthermore, two 
such forms cannot intersect each other in spacetime. The 
origin of electromagnetic radiation is explicitly traced 
back to the loss of simultaneity between the Dirac wave- 
form and its Coulomb field. 

We conclude in Sec lVIII with a short list of the existing 
data and experiments that are in line with or can serve 
as the tests for our predictions. 

The results of this work, if looked at as a launch-pad 
for further investigations, are striking in their anticipated 
mathematical complexity and physical transparency. It 
seems, however, that the former is the inevitable toll for 
the latter. The nonlinearity of the Dirac equation makes 
finding its explicit solutions a formidable task. But this 
nonlinearity is not artificial - no ad hoc nonlinear terms 
were added to the basic Dirac Lagrangian in order to 
simulate any experimentally found patterns of matter be- 
havior, symmetry, etc. On the contrary, the discovered 
generic structure corresponds to the perfectly understood 
phenomenon of localization, which is due to the local 
time slowdown, and then the loss of certain elements of 
spatial symmetry due to localization. Despite being gen- 
uinely nonlinear, these phenomena are so natural for any 
kind of wave propagation that only a minimal amount 
of information about the physical nature of the waves is 
needed to not only understand the whole picture qual- 
itatively, but even to make semi-quantitative estimates. 
In the text we also outline how the existence of the pion 
field or how the known properties of the neutrino can be 
inferred from the concept of a localized Dirac waveform. 



II. DIRAC FIELD AND RIEMANNIAN 
GEOMETRY. 

The first attempts to bring the Dirac equation into the 
framework of General Relativity (GR) was made by V. 
Fock @ and H. Weyl [1] in a series of papers in 1929. This 
study (and many other studies of that year) was in line 
with the basic concept of Einstein's GR that, in the local 
limit (inertial reference frame) , one has to reproduce the 
results of special relativity; it was established earlier that 



3 



spinors do indeed provide a linear representation of the 
Lorentz group. Somewhat later, E. Cartan pointed to a 
insurmountable difficulty - there are no representations 
of the general linear group of transformations GL(4) that 
are similar to spinor representations of the Lorentz group 
of rotations. Cartan stated the following theorem , which 
vetoed spinors in Riemannian geometry: 

" With the geometric sense given to the word "spinor" it 
is impossible to introduce spinors into classical Rieman- 
nian technique; i.e., having chosen an arbitrary system 
of co-ordinates x^ for space, it is impossible to repre- 
sent spinor by any finite number of components ipi such 
that ipi have covariant derivatives of the form ipi-,^ = 
dfiipi+^i u ipj, where P? are determinate functions of x 11 . 
"1 

Of these two underscored reservations of Cartan, the 
first one was investigated by Ne'eman et al 6], who pro- 
posed to overcome the veto by resorting to the infinite- 
dimensional representations of the Lorentz group. The 
present study explores the window, which is left open by 
the second reservation. As long as natural coordinates 
for the Dirac field are anholonomic (also in the sense 
of the theory of dynamical systems) the "connections" 
P] and, eventually, the metric tensor g l _ LV {x) appear to 
be functions of the Dirac field (which are determined by 
Eqs. glU, P~T!))) and g^J), (TC^l below) and 

not determinate functions of x^ . 

In this section we set the stage by demonstrating that 
the two key issues of geometry, direction and distance, 
can be separated in a "physical way" by associating the 
field of directions with the Dirac field of matter. The Rie- 
mannian metric will then be associated with the propa- 
gation of signals. 



A. Algebraic properties of the Dirac Field. 

All observables associated with the Dirac field are bi- 
linear forms built with the aid of Dirac matrices a 1 and (3, 
which satisfy the commutation relations ( a a ~ (l,a*); 
a = 0,1,2,3; i= 1,2,3) 

a a (3a b + a b f3a a = 2/3r 1 ab , (2.1) 

where rf h — diag(l, —1, —1, —1). We begin with a review 
of the properties of the Dirac field ip(x) which hold at 
a point, without a precise definition of the coordinates 
x. For now, ip will stand for a column of four complex 
numbers ip a . 

There are sixteen linearly independent 4x4 Hermitian 
matrices all of which can be constructed from the four 
matrices a 1 and (3. The Dirac matrices, pi { pi — (3, 
ps = —iQ V a 2 o? , p2 = —if3p3), and Oi — p^a 1 satisfy the 
same commutation relations as the Pauli matrices, and 
all a matrices commute with the p matrices: CiOfc = 
Sik + itikio-k, PaPb = Sab + itabcPc o~ip a - p a o~i = 0. The 
matrices — p% and p\ are commonly known as 7 s and 



7 , respectively . Below, these matrices are used in the 
spinor representation, 

° i = (0 -Ti)' Pl= (l Oj' 

where 73 are the 2x2 Pauli matrices. Employing the 
Dirac matrices, we can define the four components of the 
"vector current", j a = ip + a ar ip = ipj a ip, the four compo- 
nents of the "axial current", J a = ip + p^a a ip = ip-f 5 -f a ip, 
the "scalar" S = %p + p\tp = ifnj) and "pseudoscalar" V = 
ip + P2ip = — «'07 5 V'i ancl the six components of the skew- 
symmetric "tensor" M ab = {i/2)ip + [a a p\a b — a b pia a ]ip. 
The similarity of these quantities to the Lorentz ten- 
sors can be verified in a purely algebraic way. Indeed, 
if the Dirac field ip is transformed by means of a sub- 
stitution ip — > Sip (or the matrices are transformed as 
a" — > S + a a S, etc.), with the matrix S depending on 
four complex parameters, 

fx y x _ (<* 0\ 

s -{o (A+r 1 ;- x -[ 1 s)' 

det|A| = aS - f3~f = 1 , (2.2) 

then the components of j a , J a , S, V and A4 ab experi- 
ence a four-dimensional Lorentz rotation at angles which 
are uniquely determined by these parameters. For exam- 
ple, if we take a = e"^/ 2 , (3 = 7 = 0, and S = e*^ 2 , 
then the transformation S = e~ l ^ a3 ^ 2 is unitary and the 
components of j' a — ip + S + a a Sip are 

j'°=J°, j n =j 1 cos<l>-j 2 wx<l>, 
i /2 =i 1 sin^ + j 2 cos0, f 3 =f, (2.3) 

which corresponds to the rotation of the vector j a at an 
angle <p around the axis "3" . In exactly the same way, 
if we take a = e~ v ^ 2 , (3 = 7 = 0, and 6 = e v ^ 2 then 
S = e ~ va3 / 2 ; the components of j' a will be 

j'° = j° cosh j) - f sinh 77, j' 1 = j 1 , 

f = f, f = -i° sinh 77 + f cosh 7?. (2.4) 

This transformation of ip corresponds to a Lorentz boost 
in the "third" direction and is not unitary. Similar cor- 
rect relations are immediately verified for the scalars iS' 
and V' ', the vector J' a — ip + S + p^a a Sip , etc. Therefore, 
for example, the quantities 

f = Vabff = ft - ?, J 2 = VabJ a j\ 

j-J = Vabfj" = Jo Jo - JJ (2.5) 



We consciously refrain from using the anti-hermitian matrices 
7 1 = pio 1 and the Pauli-conjugated spinors ip = ip+pi. In their 
terms, the formulae of parallel transport (see Appendix A) would 
be much less transparent and unnecessarily complicated. 
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are invariants of the tp — > Sip transformations. Notably, 
the Minkowski signature matrix of Eq. (|2.1[) . n ab = r)u)8%, 
and its inverse n a b, n ab rjb c = ^ , came up here in a purely 
algebraic way, without even mentioning the Lorentz sym- 
metry and we will use it right away in order to preserve 
the usual convention about contraction of repeated up- 
per and lower indices. It is also just an algebraic exercise 
to check that K 2 = j a j a = -J a J a = S 2 + V 2 > 
and that j ■ J — 0. The latter relation means that if 
the vector current of the transformed field is of the form 
j a = (1Z, 0) then the axial current can only be of the 
form J a = (0, J) and that it can be further "rotated" 
to J a = (0,0,0,72.). Therefore, at a generic point x, the 
vectors e^(x) = j a /lZ and efy(x) — J a /1Z are the or- 
thogonal timelike and spacelike unit vectors, respectively, 
and they can be reduced to e? N = Sq and e%s = 5§. (In 
what follows, the symbol (=) is used in equations that 
imply such a particular reduction.) 

* ab 

The components of the tensor A4 and its dual M = 
(l/2)e abcd M cd are 

M m = K t = <i/j + P2 a^, M = e 0ijm K m 

M H =Li = tP+pKX^, M ij = e°^ n L m . (2.6) 

Because M ab Mbc = {L ■ K)5®, these two tensors can be 
used to build two couples of vectors which are spacelike, 
orthogonal to and e? 3 j and to each other, 

E c = { 3 a /K)M ab [5 b + J b J c /K 2 ] = (0, K U K 2 , 0), 

E c = (J a /K)M ab [S b - j b ] c /n 2 ] = (0, K 2 , -K u 0), 

H c = (J a /K)M ab [S b c - j b jc/K 2 ] = (0, -La, L u 0), 

H c = {f/U)M ab [S b + J b J c /K 2 ] = (0, L u L 2 , 0). (2.7) 

A full set of easily verifiable identities between invari- 
ants of the transformations (|2.2p is given by 

n 2 = j a f = -j a j a = s 2 + v 2 , j a f = o, 

S 2 ~V 2 = L 2 -K 2 , SV = L-K. (2.8) 

The scalars allow for the following parameterizations, 

S = 72cosT, P = 72sinT, (2.9) 

where both 1Z and T are functions of the Dirac field. It 
is important that the absolute values of S and V do not 
exceed 1Z. A similar observation is true for the second 
line of Eq.dH]), 

S 2 - V 2 = L 2 - K 2 = V? cos 2T, 

2S7' = 2L-A ? = 7e 2 sin2T. (2.10) 

Concluding the discussion of the algebraic properties 
of bilinear forms of the Dirac field at a generic point x, let 
us introduce, along with the orthogonal system e^J^], a 

reciprocal (in algebraic sense) system e£ [tp] 

a a 



where we assumed that det|e^| ^ 0. Then, a simple 
algebra verifies that the objects 

g a p = V 0b e^ef, 9 aff = V ab ef a) e% (2.12) 

can be used to move the Greek indices up and down, for 
example, 

_ (a) _ rn _ 

9af3 e (hj ~ e (a)a e p e (f,) — b e (a)a — e a(b)- 

It is also evident that the repeated upper and lower Greek 
indices are contracted. 



B. Dirac currents and Riemannian geometry. 

From now on, we look at the ip a {x) as the physical 
Dirac field, the continuous functions of the arbitrarily pa- 
rameterized points = [x 3 ) of the spacetime. 
So far, we have verified that the algebraic structure of 
bilinear forms of the Dirac field naturally contains an or- 
thogonal quadruple of unit Lorentz vectors at a generic 
point, thus defining spacetime directions at that point. 
In a sense, linear transformations (|2.2[) of the Dirac field 
generate a group of homogeneous linear transformations 
for vectors, thus associating with every point a local cen- 
tered affine space. Vectors of this quadruple are thought 
of as smooth functions of ip(x) and it is tempting to im- 
mediately treat it as a quadruple of the vector fields. 
But these transformations of the Dirac field have noth- 
ing to do with the general transformation of coordinates, 
which are arguments of ip{x). For a given fixed A, we 
can consider x x = const as the equation of a coordinate 
hypersurface and the lines along which all coordinates 
but x x are constant as coordinate lines. Tangent vectors 
of these lines (which are gradients of the linear function 
ip(x) — x x ) are = dx^/dx x = S^y, therefore, this co- 
ordinate system is an holonomic one, but it has no metric 
and there is no way to determine if its coordinate lines 
are orthogonal. One may replace x^ by smooth functions 
of other coordinates y M , x^ = /^(y), thus redefining co- 
ordinate lines and surfaces, but such a change does not 
alter ip(x{y)) and has nothing to do with affine Lorentz 
transformations (|2.2[) . To bridge the gap between the ab- 
stract field of directions determined by the Dirac field and 
the given above definition of the holonomic coordinates, 
it is necessary to know in advance that four systems of 
differential equations (for the unknown 

da; dx 1 dx 2 dx 3 

for congruences of lines (labeled by the ordinal numbers 
(a) ) are solvable and thus determine a coordinate net. 
In other words, if these equations are integrable, then the 
system 

dx a = e? a) ds a , ^ = 0,1,2,3, (2.14) 
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will represent lines, which at every point x have a deter- 
minate direction e" a j( , 0), and only one line of the con- 
gruence (a) passes though each point in spacetime. The 



tetrad < 



(a) 



will be a Lorentz vector and a coordinate vec- 



tor. A change of the vector variables x^ = f^(y) in 
Eq. (|2.14[) will result in the transformation of the differ- 
ential dx^ , 

Br a 8r a dii a 

dX ~ dy V ' dy 8x0 - a ' 

and the equation for the same congruence in new coordi- 
nates will read as 

dy" = el a) {y)ds a , efa(y) ^ |g • ef o) [# (2.15) 

A new element here is that tangent vectors depend 
on the coordinates via the coordinate dependence of the 
Dirac field which, in its turn, is constrained by the equa- 
tions of motion. Therefore, the problem of integrability 
of Eqs. ()2.13|l - ()2.15|) cannot be addressed solely within 
Riemannian geometry; at least some properties of con- 
gruences must be controlled by the Dirac equation. For 
"temporal" and "radial" congruences, the Dirac equation 
indeed yields a set of constraints with a clear physical 
meaning. The properties of congruences of angular arcs 
(including their symmetry), in general, not only explic- 
itly depend on particular solutions of the Dirac equation 
but there may even be no meaningful holonomic coor- 
dinates associated with these arcs. Nevertheless, even 
keeping such a difficult perspective in mind, let us con- 
sider all four tetrad vectors as contravariant vectors of 
Riemannian geometry. 

Traditional approaches assume solving the Dirac equa- 
tion in a determinate metric field g^vix) of spacetime. 
The polarization properties of the Dirac field prompt the 
opposite direction of thinking. Namely, the field ip(x) 
must be the solution of the Dirac equation, which ex- 
plicitly depends on a resulting metric g^[ip(x)] given by 
Eqs. (|2.12[) . In such a context, the Minkowski form of the 
metric in the local limit is associated not with an imag- 
inable local inertial frame but rather with the algebraic 
properties of the Dirac field and (complementary to the 
latter) the hyperbolic character of the Dirac equation. 
Thus, it is possible to overcome Cartan's veto in two ma- 
jor points. First, there is no arbitrary coordinates for 
spacetime (modulo a trivial change of variables). Sec- 
ond, the connections, T, are no longer determinate func- 
tions of x; they become functions of the Dirac field. In 
this framework, as is shown below, the hypersurfaces of 
a constant temporal coordinate naturally emerge; their 
existence is a prerequisite for the quantization of the on- 
mass-shell Dirac field in curved spacetime. The proper 
time slows down in domains of a higher matter density, 
which points to self-localization as an intrinsic property 
of the Dirac field. This effect also clarifies the nature of 
electric charge and of charge asymmetry of the empiri- 
cally known stable matter. Along with localized matter, 



there always exists a preferred system of orthogonal con- 
gruences determined by the internal polarization struc- 
ture of physical objects. In general, it can be considered 
as the net of the anholonomic coordinate system Q. 

Let us follow the key idea of intrinsic geometry to asso- 
ciate tensor fields with mutual invariants of tensors and 
parameters (tangent vectors e^x(af)) of a system of con- 
gruences. Furthermore, let us read Eqs. (|2.11[) and (|2.12p 
as 



r, ab e^(x)e^(x)^g^(x) 1 



(2.16) 



and consider this gufi(x), as a primary choice of the space- 
time metric. Only by virtue of Eqs. (|2.13|) can we trans- 
late the first of equations Eq. (|2.16p into 

ds 2 = 7]abds a ds b = g l/tl (x)e^ a) (x)e" b) (x)ds a ds b 

= g v „(x)dx^dx u , (2.17) 



and reconcile Eqs. (|2.12p (inspired by the algebra of the 
Dirac matrices) with the measure of length postulated in 
Riemannian geometry. 

When e£i is a vector with the law of transformation 
(|2.15[) and g v ^(x) is a tensor (not necessarily determin- 
ing a metric) then the covariant derivative V^e!^ with 
respect to g v ^ is also a tensor [8|. Therefore, one can 
introduce a system of invariants (the Ricci coefficients of 
rotation of a system of congruences) 



(2.18) 



For a given (c), six parameters Lu a b c ds determine an in- 
finitesimal rotation of the pyramid of tetrad vectors in 
the "plane" (ab) when the vertex of the pyramid is dis- 
placed by ds along a line of congruence (c). Equation 



(2.19) 



is the inverse of lpT8]) . Using Eqs.jHHl) and |2T9|) . it 
is straightforward to check that if g v ^{x) has the form 
(|2.16p then V \g vil = 0. Consequently, the vector con- 
nections r^„ coincide with the Christoffel symbols of the 
metric g ufl . 

In an ideal geometric world (i.e. when all four holo- 
nomic coordinates exist) the necessary conditions for in- 
tegrability of Eqs. (|2.19p are given by the Ricci identities 



(V M V A - V A V A 



(2.20) 



where R av ^\ is the Riemann curvature tensor. These 
equations can be cast in the form, 



a v ^ A o d 

-(a) e {b) e (c) e (d) tt <?vti>> — U abcd, 



(2.21) 



where 



Rabcd = ddOJabc ~ 9 c UJ a bd (2.22) 
Vf [Ufad^fbc — LOf ac LOfbd + ^>abf(^>fcd ~ ^fdc)), 

f 
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is a system of invariants, which is then known as the 
tetrad representation of the Riemann tensor. Since at 
least some of the Ricci coefficients of rotation will ap- 
pear to be functions of the Dirac field, this dependence 
will be carried through onto the Riemann and Ricci ten- 
sors. The Einstein equations for the metric field g V p,{x) 
that describes motion of macroscopic objects may appear 
to be descendants of the constraints stemming from the 
Dirac equation. 

To summarize, if in spacetime, with arbitrarily chosen 
holonomic coordinates, x 11 , the Dirac field ip(x) is de- 
fined and at each point the 16 quantities, e^[ip(x)], are 
computed (along with the algebraically reciprocal system 

eS [^0*0]) then the metric g^^(x) of spacetime is given 
by Eq. pTIo]) and the interval by Eq. ([2T7|) . This metric 
depends on the Dirac field and is not defined a priori. 
From the physical perspective, its existence seems to be 
a privilege of exceptional solutions rather than a rule. 

It is important to realize that the material Dirac field 
defines a system of the unit vector fields {x) - there- 
fore, the effect of such a matter- induced metric should be 
equivalent to a long-range interaction between localized 
objects. 



III. DIFFERENTIAL IDENTITIES FOR 
TENSORS. 

In order to find limitations on the metric of spacetime, 
which can host the localized configurations of the Dirac 
field, we begin with the examination of various identities 
that are consequences of the Dirac equation. The ques- 
tion is, whether differentials of various bilinear forms of 
the Dirac field, which are considered as the physical ob- 
servables, can be translated into covariant derivatives of 
tensors. We use this question as a test of the roots of the 
discrepancies that could have led to Cartan's veto. It ap- 
pears that these discrepancies correspond to the clearly 
understood physical processes. 

Following Fock and Weyl, we postulate that the equa- 
tion of motion of the Dirac field and its conjugate are 

a a D a ip — —impiip, (3-1) 
?/> +< D^a a = irmp + p u (3.2) 

where the covariant derivative D a ip = (d a — T a )ip of 
the Dirac field is defined in Appendix [XJ The object 
D^ip = e^D a ifj = (dp — Tfj^tp will be used only as a sym- 
bol, since it has no clear geometrical meaning. The mass 
parameter in these equation is a priori arbitrary. Be- 
cause the Dirac field has the property of self-localization, 
every stable localized waveform will determine the corre- 
sponding value of to. 



A. Identities for vector and axial currents. 

From the equations of motion (|3.1[) and (|3 . 2[) one im- 
mediately derives two well-known identities. One of 
them, 

D a f = V = J=d? lV=g-i>+a^} = 0, (3.3) 

clearly indicates conservation of the vector (probability) 
current of the on-mass-shell Dirac field, while the second 
one indicates that the axial current is not conserved, 

D a J a = V = 2mV, (3.4) 

and has the pseudoscalar density as a source. The sig- 
nificance of Eq. (|3.4p is due to the pseudoscalar density V 
on the r.h.s. Since V is localized not less than 1Z and the 
vector J^ L is spacelike, the unit axial vector e^L defines 
the outward radial direction. The existence of such a 
direction is a distinct characteristic of a localized object. 

B. Flux of momenta: not tensors. 

Consider now a more complicated object T b a = 
iip + a a Db'il>, the Hermitian part of which is normally 
regarded as the energy-momentum tensor of the Dirac 
field 2 . Its components are interpreted as the flux of com- 
ponents iDi, of the momentum in the direction of con- 
gruence of lines of the vector current j a . Because the 
vector current is timelike, this tensor is well-suited to 
describe the flux of momenta carried by massive parti- 
cles. The Lagrangian Lo of the Dirac field vanishes for 
the on-mass-shell configurations and T 6 a does not have a 
diagonal term, —Ld8%, which could have been responsi- 
ble for the flux of momentum in the spacclikc direction 
(e.g., the pressure). Since the spacelike radial direction is 
controlled by the axial current, we are led to consider an- 
other object, the stress tensor P 6 a = iip + pz^Dbifj, which 
accounts for the flux of momenta in the radial direction. 
For stable localized wave forms, there must be no flux of 
any observables in the spacelike outward direction. How- 
ever, if we decide to investigate a particle's Lorentz con- 
traction as a dynamic process or the decay of a long-lived 
waveform (considered as a particle), then we are led to 
consider the spacelike flux of momenta due to the "phase 
shifts" inside the wave form. Regardless of how adequate 
this intuitive physical interpretation of T£ or P b a is, or 
even without any physical interpretation, they both can 



2 The reader should not be confused by how the standard name 
"energy-momentum tensor" is used. In the context of the present 
work, the invariants T£ and P£ are the auxiliary objects. We are 
interested only in identities that can be derived from the Dirac 
equation in tetrad form and then translated, if possible, into the 
tensor form. Only Hermitian part of these objects enters the 
equations that allow for a physical interpretation. 
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be used to derive various useful identities , which allow 
one to compute the rotation coefficients aj a b c as functions 
of the Dirac field and thus constrain the possible metric 



(I2.17p . In this section, we study T b a in detail. The stress 
tensor P b a is studied in Appendix B along the same guide- 
lines. 



I 

One would expect that the absolute differential of T ab , being computed according to the Leibnitz rule, will be as 
follows, 

D c T a b = d c T a b — UJ a dcTdb — UbdcTad = V ' c T ab ■ (3.5) 

If this expectation turns out to be justified then the usual covariant derivative will be immediately reproduced as 

d x T^- K x T Ufl - r£ A T w = e^e* V c T a6 = V A 7V (3.6) 
Contrary to the expectation of (|3.5p . the answer reads 

D c [^ + a a ~S b ^] = d c [^ + a a ^b^} - V + [r+« a + a a T c b Tp = d c [^+a Q ^] - uj adc ^ + a d ^ b ^, (3.7) 

with the last term of Eq. (|3.5[) missing, and no hope to recover the full expression (|3.6[) of the covariant derivative of 
the tensor! Contracting indices a and c we arrive at the expression, 



D a [4> + a a ^ b ip} = d a [^ + a a ^b^]+uj acc ^ + a a ^b^ = 



-g dx v 



(3.8) 



which is exactly what one may wish to have as the l.h.s. of a conservation law for the energy-momentum of the Dirac 
field . The missing term is exactly the one that does not let one interpret equations like Vo-T^ = as conservation of 
anything. However, at the moment, a covariance can not yet be explicitly visible; it may occur that the r.h.s. of an 
expected conservation law recovers the covariance of the resulting identity as a whole. The r.h.s. must be determined 
using the equations of motion. Let us first rewrite the l.h.s. of (|3.8p as 



By virtue of the equations of motion (and due to the 
Leibnitz rule) the last three terms exactly cancel out and 
the final result is 

D a [i4,+a a ^ b ip] = i^ + a a a ~S b - ^b^aW- (3.9) 

Using Eqs. (|A.10[) and (|A.11[) to separate the terms with 
and without derivatives in the commutator and compar- 
ing with (|3.8[) we find that the Dirac equation yields the 
identity: 

d a T ab - oJcacTab = UbcaTac- u cab T ac - a a TD ab ip , 

(3.10) 

where the commutator Vi ab does not contain derivatives 
of ip. After moving the term ujbcaT ac from the right to 
the left, the l.h.s. becomes, according to (|3 . 5|) and (|3.6[) . 
the system of mutual invariants of a usual covariant di- 
vergence of the tensor T° and congruences e^. It can 
be transformed either into coordinate form (|3.6[) . which 
has no explicit dependence on tetrad vectors or into a 
non-coordinate form. Unfortunately, this coordinate in- 
dependence of a fragment of identity (|3.10p is useless, 
because there remains an abnormal term iO cab T ac on the 
right. Being translated into a coordinate form, this term 



becomes (y \e v a )e%T^ l . It explicitly depends on tetrad 
vectors (on how the coordinate lines are bending). 

The abnormal term Lo ca bT ac in Eq. (|3.10p enters another 
identity that follows from the Dirac equation. It arises 
after contracting indices a and b in Eq. (|3.7p . 

D c [i; + a a ^ a ijj}=d c [ijj + a a ^ a ij}~ Lo abc ^ + a b ^ a ^. (3.11) 

Eq. (|3.1ip reveals one more inconsistency, which is similar 
to the one observed in Eq. (|3.8p . The quantity D c T a a is 
derivative of the trace of a tensor, i.e. of a scalar. How- 
ever the result has an additional term with a connection, 
which is one more piece of evidence that the quantities, 
T a b , are not the invariants of a tensor. By the same to- 
ken, the l.h.s. of Eq. (|3.1ip is not a covariant derivative 
of a scalar. 

By virtue of the Dirac equation, the first term on the 
r.h.s. of (|3.1ip becomes d c [—im4> + piip}- Alternatively, 
one can immediately use the equations of motion on the 
l.h.s. and only then differentiate, 

D c [4' + ^ a ^] = -imD c [^ + Pl ^} 
= -imd c [ip + p!i/j] + imip + [r+ Pl + piT c ]tp . (3.12) 
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Comparing the last two equations and using (|A.4[) we 
finally get 

toacb ■ T ca = 2mgVK b . (3.13) 

Using Eq. (|3.13[) , one can then rewrite Eq. (|3.10[) in a for- 
mally covariant form, 

V a T a v = i^+a^D^ D v ]i> + 2mgVK , (3.14) 

where = N Q and the commutator D ^D V — D V D ^ = 
[D^Dv] = — e°e(j,D a f, on the r.h.s. has no derivatives. 
For the sake of completeness we mention that the imag- 
inary part of T a v is a tensor; it is the covariant deriva- 
tive (i/2)V v j a . Because the vector current is conserved, 
the imaginary part of Eq. (|3.14p is just an identity [H, 
«V CT (V I ,j <T ) = iRov3° ', where is the Ricci curvature 
(contracted Riemann tensor of curvature). 

An attempt to make N„ = leads to the main re- 
sult of Fock's paper Q , which was derived entirely in the 
coordinate representation (using as a well-defined op- 
erator) and interpreted, with the reference to the corre- 
spondence principle, as the equation of a geodesic line. 
Since Eqs. (|3.9p and (|3.13p are nothing but two identities 
that follow from the Dirac equation and Eq. (|3 . 14[) is their 
sum, there obviously is a way to derive Eq. (|3.14p in one 
step, which was done by Fock (with a subtle inaccuracy). 
The possibility to set H a = can be viewed as evidence 
that the abnormal term w ca f,T QC is zero. Such an impres- 
sion is not correct. This would require that oJ ca b = 0. If 
the second term in the r.h.s. of Eq. (|3.10p is zero so is the 
first one. There remains nothing to move to the left in 
order to compile the covariant derivative of the tensor. 

Consider another possibility that the Riemannian 
spacetime, which is hosting a Dirac field configuration 
(like proton) admits an orthogonal system of coordinate 
hypersurfaces. Then the Ricci coefficients with all differ- 
ent ordinal numbers vanish, i.e. uj ca b — Wcba = 0, and the 
first two terms in the r.h.s. of identity (I3.10P just cancel 
each other. Once again, the possibility to compile the 
covariant derivative of Tab as a part of the conservation 
law is lost. We are led to conclusion that the metric of 
spacetime, which is hosting the Dirac field (and thus is 
determined by this field) does not allow for a system of 
orthogonal coordinate surfaces. 

If we formally translate the remaining terms into the 
coordinate representation, we get, instead of (|3.10p and 
()3.13p . two equations, 

= iyj+ a v[D^D„-D„D^, (3.15) 
(V^KT- = YM - KJTI = 2mgV%, (3.16) 

both carrying an explicit dependence on the tetrad vec- 
tors. The sum of these equations reads as, 

d T a + T a T v - T v T a - id e a - d e a )e"T' J 



where this dependence is apparently hidden because we 
assumed that all congruences are normal and d a e a ^ — 

^M e a = Sai^abc — ^acb)^n = 0. This equation indeed co- 
incides with (|3.14p . but only when the connection, r£„, 
is symmetric, which was not an a priori requirement. 
Since, in general, the Ricci coefficients are not zero and 
the "tensor" T afl is not symmetric (except for a plane- 
wave solution) , we cannot argue that the r.h.s. of (|3. 13|) 
must be zero for whatever reason. If, in addition to (jA.ip , 
we unconditionally required that 8S = SV = 0, then ar- 
riving at (|3 . 1 3[) we would generate controversy. 

An ad hoc choice of an orthogonal coordinate sys- 
tem (where uj a bc — Wacb) can serve only as a crude 
approximation. To be consistent, we have to replace 
UacbTca by u! a b c T ca in Eqs. (|3.10p and Eq. (|3.13p simulta- 
neously. Then, the latter can be rewritten as uobca • T ca = 
—2mgV&b, so that the symmetric (and only the symmet- 
ric) part of Tb a matters. This transmutation indicates 
that we implicitly employed the approximation of a ma- 
terial point when internal deformations, that are bringing 
an object into a new state of motion, are discarded. In 
this case, the unit vector plays the role of the 4- 
velocity of a small object as a whole. Since the first 
term in brackets in Eq. (|3.16p can be dropped, we may 
write 

^ [V=5 Re(Tt)] = e^—gfF^ v , (3.17) 

Re(T* ) = -2mgV% , (3.18) 

which is a perfect expression for the energy-momentum 
conservation complemented by the constraint (|3.16p . If 
the equation W a T a v = is considered as a prototype for 
the equation of a geodesic line (as it was conjectured in 
0) then the term T^T^ = Y v ajl T a v in it is connected 
with H M through Eq. (|3.18p . Depending on the nature of 
the physical process, this term is responsible either for 
the gravitational force or for the force of inertia. These 
forces are real and one cannot set or V to zero without 
losing them. An estimate of the coordinate dependence of 
H M yields Newton's approximation for the metric tensor. 
At large distances, we have H oc 1/r 2 , as it follows from 
Eq. (|4.18p . A startling connection of the field V with the 
localization of the Dirac field and origin of its mass is 
discussed in SeclIVI . 

The physical origin of these forces can be understood 
from another perspective. The r.h.s. of Eq. (|3.16p can be 
rewritten as m{d a S — D a S). The first term accounts only 
for propagation of the wave form considered as an object. 
The second term also accounts for the internal polariza- 
tion of the wave field. The difference between them is a 
force, which is due to internal polarization. This fact mo- 
tivates the view on coordinates, as descendants of the po- 
larization structure of the Dirac field, which was proposed 
in Sec lII Al Its dynamics are described by Eq. (|B.10p . An 
immediate consequence of the existence of an internal dy- 
namic in a localized Dirac waveform is a view of pions 
as one of polarizations of the Dirac field, which can be 
resolved, e.g., as a 27-resonance (see Appendix [Bjl . 
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IV. DIRAC FIELD AND CONGRUENCES 
OF CURVES. 

In this section, we closely follow the ideas of the in- 
trinsic geometry of Ricci and Levi-Civita as they are pre- 
sented in the monograph [9i] ; the metric properties of the 
spacetime are expressed in terms of rotations of the local 
coordinate pyramid. The main subject of the analysis 
are Eqs. (|3.13p and (|B.4[) . which are the differential iden- 
tities that follow from the Dirac equations. Eq. (|3.13[) is 
trivially satisfied only for plane waves, i.e., when V = 
and the tensor T a {, is symmetric. These solutions are 
employed in scattering theory and they do not repre- 
sent particles. In such a context, equations like (|3 . 13[) 
and (|B.4j) cannot even be derived. Unlike the commonly 
known identities (O and (EH, Eqs. ([3~T3|) and |B~4|) are 
not covariant in the sense that they explicitly depend on 
congruences, which are the physical characteristics of the 
Dirac field. It appears that these identities impose im- 
portant limitations on the properties of the metric, which 
is compatible with the localized solutions of the Dirac 
equations. These limitations are studied in the following 
section. 



A. Vector current and timelike congruence. 

The Ricci coefficients are real-valued and skew- 
symmetric in the first two indices. The tensor T a b is 
neither real nor symmetric. The r.h.s. of Eq. (|3.13[) is 
real. Therefore, the imaginary part of Eq. (|3.13p is just 
Im(T ac - T ca ) = D c (il; + aaip) - D a (ip + a c ip) = V c j Q - 
VqJc = 0, and it should be considered together with the 
equation V a j a = of the vector current conservation. 
Since \7 a jb are the invariants of a true tensor, V^j„i we 
have two tensor equations, 

V^V - V„J P = . (4.1) 

and 

V^f = . (4.2) 

The vector field j M (x) is globally timelike; its tangent unit 
vector is e*f Q Jx), j M = 1Z e? Q s, where 1Z = \J~f- > is the 
invariant density of the Dirac (spinor) matter. Therefore, 
Eq.(|4.ip becomes 

V^eS 0) - V„e<°> + Inll - e^d v lull = 0. (4.3) 

Contracting this equation with e^e^, a, b = 1, 2, 3 and 
recalling Eqs. (12.18p we find that 

UQab ~ UQba = , a, 6=1,2,3 (4.4) 

which is a necessary and sufficient condition for the con- 
gruence e^jj to be normal [1, Namely, there exists 

such a function, T(x), that the vector field \x) is or- 
thogonal to the family of surfaces T{x) — const, 

d li T{x)=f{x)ef{x), (4.5) 



where f(x) is a coordinate scalar. Contracting Eq. (|4.3p 
with we get 

Inft = e^d (0) hxTl - wwoe^, (4.6) 

where <9(o) In 72. = e? sdn In 72, is the derivative in the di- 
rection of the arc so . Contraction of Eq. (|4.3p with 
yields 

(dlnK/ds a ) = ~Lu a0 o , a =1,2, 3, (4.7) 

which indicates that congruences of lines, defined by the 
system of equations (|2.13|) . dx^/dso = , must ex- 
perience permanent bending (acceleration) whenever the 
invariant density lZ(x) of the Dirac field is not uniformly 
distributed. The spatial gradient of IZ(x) cannot vanish 
for any localized state. Even more, the congruence of 
lines of the Dirac current is not a geodesic congruence, 
since, for geodesic lines, the vector of geodesic curvature 
would have vanished, i.e., u> a oo — 0. 

Additional information can be extracted from Eq. (|4.2|) . 
From definition (|2. 19[) it follows that 

V„e^ 0) = -(dlnn/ds ) = J2v(a)U aa ■ (4.8) 

a 

Hence, we can rewrite (|4.6p as 

CV In K = -4°) V(a)U aa ~ WWO^, (4.9) 
a 

which shows that the r.h.s. of Eq. (|4.9p . which contains 
only geometric objects, is a component of a gradient. 
Together with condition (|4.4p this constitutes a necessary 
and sufficient condition that the function T(x), defined 
by Eq. (|4.5p . is an harmonic function [8j, 

□ T = Q^V^JT = 0. (4.10) 

We may take the parameter t of T(x) = t = const as a 
definition of the world time. For the harmonic function, 
T(x), the conditions of integrability for the system (|4.5p 
of partial differential equations reads as @ 

d^lllf = -B^ Y V(a)^Qaa ~ ^600^ ■ 
a 

Comparing it with (14. 9p we find that f(x) — TZ, so that 
the world time t and the "proper time" Sq are related by 

dt = 1Zds = cfeo/V5oo ■ (4-11) 

Hence, we can draw the major conclusion that: The 
proper time, sq, flows more slowly than the world time, t, 
whenever localized Dirac matter is present. If Dirac mat- 
ter is in a stable configuration (on mass shell), then there 
is a well defined time and one can consistently speak of 
a (quantum) state of the Dirac field. 

Since the congruence appeared to be normal, the 
hypersurfaces T{x) = t = const represent space at dif- 
ferent times t. The three other vectors e^Ax), i = 1, 2, 3 
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of local tetrad are spacelike, orthogonal to e^(x) and 
thus belong to such hypersurfaces (by the definition, 



d a T 



0). The interval becomes as 

goadt 2 + g lk dx % dx k . 



ds z 



(4.12) 



Accordingly, eft 



M - 



(0) 



Equation ()3.3p of the vector current conservation now 
reads as 



0. 



This can be recognized as the condition for the coordinate 
x M to be harmonic, 



dx^ 



0. 



This can be read in different ways. First and foremost, 
it expresses the curvature of the two-dimensional sur- 
face (si,S2) of angular coordinates via the local param- 
eter sinT = V /1Z of the Dirac field. Vice versa, once 
geodesic curvatures W131 and W232 are known in advance 
(e.g., from an alleged symmetry, experiment, etc.) then 
T(ir) is known as an explicit function of spacetime coor- 
dinates and there remains no freedom of "crural" trans- 
formations, like in Eq. (|A.8[) . 

Second, the l.h.s. of (|4.16[) must be a well-defined 
(4.13) geometric object (at least when the congruence e^ is 
normal and the radial coordinate is well-defined). In 
this case, we must have D a T = because the covari- 
ant differential operator D a is defined only by its action 
on the Dirac field. Consequently, by virtue of Eq. (|A.8|) . 
D a T[ijj] — d a T + 2gK , we have 



which is specified for the normal coordinate ip = T = x 
(see, e.g. [Tg], §41). From P~HS|) , there follows one more 
(very intuitive) form of the current conservation, 



(4.14) 



where g = det|gjfc| is the determinant of the spatial met- 
ric. This form means that when the density 1Z grows and 
local time slows down, then the measure of space volume 
shrinks. Since the variation of 1Z and time slowdown 
both are intimately connected with acceleration, the last 
equation unites them and the Lorentz contraction of a 
localized object in one physical process. One should not 
even refer to a spacelike interval between two events on 
the opposite sides of an elementary object. 



B. Axial current and spatial part of metric. 

The axial current J 11 is spacelike and orthogonal to the 
vector j 1 *. According to Eq. (|3.4|) . the axial current has a 
source 2rriP, which is localized together with the invari- 
ant density 7Z. Since there is no flux of vector current 
in this direction (the amount of matter inside a closed 
surface remains the same), we associate the radial direc- 
tion e? 3 \(x) with the axial current, J 11 — TZeLy Then 



Eq-dSJ takes form 



V„e 



C & (3) ^ °(3)^ 

On the one hand, by definition, 



d^lnTZ = 2mV/TZ = 2m sin T . (4.15) 



M c (3) 



, V(a)^>3aa — ^>300 ~ ^311 ~ ^322- 



On the other hand, according to Eq. (|4.7p , we have 



'(3)^ 



d^lnTl = d\nll/ds 3 = -w 30 o- 



Substituting these expressions into Eq. (|4.15p we obtain 
an extremely important relation, 



^131 + ^232 = 2m sin T 



(4.16) 



2gK = -d a T, 



(4.17) 



i.e. the field N a must be a gradient. If the congruence 
e^ 3 % is not normal, then any symmetry of any explicit 
solutions of the Dirac equation with respect to arcs si and 
S2 should be considered a dynamical internal symmetry. 

Third, for a concave surface the curvature is positive so 
that < T < 7r/2. For the normal orthogonal spherical 
coordinates we have W131 = W232 = 1/r and if such a 
coordinate system were possible we would immediately 
know that 



T[ip] — arcsin(l/mr), mr > 1 
2<?N 3 [V] = -d r T = ' 



1 



(4.18) 



Obviously, this simple formula cannot be exact; rather 
it predicts the correct asymptotic behavior at large dis- 
tances. 

Fourth, the condition sinT(x) < 1 defines the mass 
parameter m as the upper limit of a possible curvature, 
which is, in fact, a definition of mass from the perspective 
of the internal structure of a Dirac particle. (In spherical 
case we would have mr > 1; the radius must exceed the 
Compton length!) This result is also in agreement with 
the kinematic Lorentz contraction of special relativity. 
An accelerated particle is Lorentz contracted and both 
7Z and the maximal curvature become oc — v 2 . 

In order to facilitate further analysis of the real part 
of Eo. (|3.13p , let us rewrite it's l.h.s in terms of the axial 
current. Let us use the dual representation of the axial 
current as e stua J a = iip + a s pia t pia u tp, (s,t,u,^), and 
differentiate it. The result reads as 



7-) stua rj 
da 



Du^+aViaVia"^)- (4.19) 



If we extend here the sum over all values of u (this sum 
vanishes by virtue of the equations of motion) and sub- 
tract the terms with u = s and u = t, the result will 
be 

D u e stua J a = -i^+a. s ~6 t i> + i^+a 1 ^^ 
-iip+D'^a*^ + iip + D~^ a s ip , 
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where the r.h.s is four times the anti-symmetric Hermi- 
tian part of the energy momentum tensor. Therefore, the 
real part of Eq. (|3.13p reads as 

(l/4H cb • e acst ■ V s J t = 2mgVK h , (4.20) 

and can be immediately recognized as the dual to 
Eq. (|B.5|) . derived for the stress tensor, 

(l/2)u) acb V c J a = -2gmSX b . (B.5) 

These two equations clearly indicate that any motion 
of the Dirac field follows the path of a helix. The un- 
derlying reason for such a motion is the non-vanishing 
curl of the spacelike axial current. The boost woi a in the 



At this point, we can conclude that the parameters Hf, in 
the connection Tb (|A.3j) are totally defined by the bend- 
ing of the system of congruences. Despite that fact that 
the Lagrangian for the Dirac equations (|3 . 1 1) and (|3.2[) 
includes the term j7 a H a , which can be interpreted as an 
interaction between the axial current and the field H a , 
it cannot be viewed as an independent field that is gov- 
erned by an additional equation of motion. (Otherwise, 
such equations must be invented, which we, so far, tried 
to avoid.) 



C. The case of the normal radial coordinate. 
Qualitative consequences of the localization. 

Even for the localized waveforms, the existence of the 
surface of a constant distance from a center is not given 
gratis. Such a surface must be orthogonal, at every point, 
to the tangent vectors e^ 3 , of the congruence of lines of 
the axial current. Unlike the previously studied case of 
the timelike congruences e^, the corresponding condi- 
tions for integrability do not universally follow from the 
equations of motion. Most likely, the radial coordinate 
cannot be normal. However, sometimes (mostly for long- 
lived particles) empirical data may hint that such a nor- 
mal hypersurface of a constant distance S3 from a cen- 
ter, which is spanned by the "angular" arcs (s\, S2) with 
tangent unit vectors (|2.7|) may be a good approximation. 
This is what we intuitively expect in a one-body problem 
and we have to verify that this assumption is consistent 
with the equations of motion and the established earlier 
constraints. In what follows, we consider Eqs. (|4.18p as 
the criterion of the spherical symmetry and try to profit 
from the fact that Eqs. (|4.2ip and (|4.22p significantly sim- 



direction Si is inevitably accompanied by the spatial ro- 
tation e l; ! k Ljjka in the plane perpendicular to s^. In plain 
words, the Dirac field cannot be accelerated without caus- 
ing a rotation thus behaving as a (relativistic) system of 
inertial navigation. 

In Eqs.(j420]) and (fB~5|) . V s J t = tf s f( t) V^Jv Since 
= T^e^-j , we further have 

D s J t = ftefoefoV „e< 3 > + Sfe^K 
= n[LU 3ts + 6?(d\nK/ds s )} . 

Finally, using Eqs. (j4.7|) and (j4.8p . which define lo s oo &nd 
^Oaa as the functionals of the Dirac field, we arrive at 



(4.21) 
(4.22) 

I 

plify under the assumption that the congruence e^, is a 
normal congruence (which, possibly, can be a first ap- 
proximation in a sequence of iterations). This will also 
allow us to qualitatively understand the trends in criti- 
cal behavior of the tetrad vectors near the limit surface, 
sin T — > 1 , and rediscover some well known properties of 
matter (which cannot be done in a picture of matter as 
plane waves). 

Since the congruence is set normal, there should 
exist a function C(x) such that 

fl/(j) = ((i)e( s )(i), (4.23) 

where is a coordinate scalar. The hypersurfaces 

C{x) = r = const are the surfaces of radius r, i.e., 

dr = (ds 3 = rfs 3 /V933 • (4.24) 

From the integrability condition for Eq. (|4.23p it is 
straightforward to derive the equations (which are similar 
to equations for the the function f(x) (cf. (|4.5p ) 

- InC = e£%) InC - oj a33 e^ , (4.25) 

(dlnC/ds a ) = co a33 , (a = 0,1, 2); (4.26) 

but, we have no constraint that would express ( as a 
function of the Dirac field. For the normal congruence 
we have oj 3a b — 0J 3 ba f° r a, 6 ^ 3, a 7^ 6, as a nec- 
essary and sufficient condition [j| Q and, consequently, 
the first term in brackets in Eqs. (|4.2ip and (|4.22p sim- 
plifies, uJ st bU 3st = -uj 03b uJo33 + ^I3b^i33 + W23&W233 and 

L °acb£ aCStuJ 3s t — ^126^033 + w 026 w 133 + w 01fc w 233- in t ne 

second term, the sum includes only a, c = 1, 2. Because 



I 

(1/4) uj acb ■ e acst [Lo 3ts - <5 3 w s00 - SfS° s ujoaa} — tti sin T ■ 2 5 H 6 , 
(1/2) Lo s tb ■ [Lo 3ts - StLUsoo - SfS°uj aa] = -mcosT ■ 2gH 6 , 
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we assume a stable object, the third term in brackets can- 
cels, dln(K(;)/ds » 0. Hence, the system of Eqs. (|4~2Tj) 
and (|4.22|1 can be cast as 

{uJiLo 02 b - w 2 cj ib) = 2mgK b sinT, 

(cjicj 3 i fc + w 2 cj32b) = 2mg^ b cos T, (4.27) 



where 



_ Vjaa + Uj33 _ 1 d\n{(/TZ) . _ 
Uj ~ 2 ~2 ds< ' J ~ ' 



Then, Eqs. (|4~2T)) with 6 = 0, 1, 2 (and H = N 1= N 2 = 0) 
yield a set of six equations, 



_ 1^320 _ ^011 _ ^012 _ W321 _ ^232 
<^2 W310 W021 W022 ^131 w 312 



±1, 



(4.28) 



where the rightmost equation immediately follows from 

Cl>3 12 = Wl31^232 = l/?" 2 and U 131 = UJ23 2 = V r - In 

spherical case, H3 is given by Eqs. (|4.18p and, conse- 
quently, 



2wiWQ23 



r 2 y/m 2 r 2 

U>2 — WQ13 = TWQ23; W312 



= , ^l(^313 ± W323) — 

1 r" 3 



W321 = ±l/r. (4.29) 



As one could expect, in the case of normal radial congru- 
ences, there is a full symmetry between congruences of 
arcs dsi and ds2- At large distances we generally have 
mr 3> 1 so that spatial rotations dominate. Vice versa, 
near the inner boundary mr — 1 the boosts ujQi^dss and 
^023^53 m tangent directions (as well as boosts Wo3i<isi 
and ^032^52 in radial directions) become infinite. When, 
starting from a generic point Xq, we move along lines 
of congruences e^(x) approaching r <~ m , then lo- 
cal tetrad rotate (with respect to the e^, (xq)) in such a 
way that all directions become nearly lightlike, so that 
the tangent velocities Vi — s, — > c. These observations 
explain the formally derived inner boundary r = 1/m 
(generally, sinT = 1) of the Dirac particle as the caustic 
of the lines of the Dirac currents. 

In fact, we have two interconnected mechanisms of the 
time slowdown (due to the amplified 1Z and because the 
vector currents tend to approach the lightlike directions), 
which cannot be separated. From the perspective of an 
"external observer", the time flow literally stops at the 
critical surface of a stable Dirac waveform. Therefore, the 
sharp interaction of the deeply inelastic scattering always 
resolves an apparently static object in a random config- 
uration determined by the foregoing causal evolution (cf. 
discussion of evolution equations of QCD in Ref . [13l] ) . 
Certain patterns of symmetry observed in such processes 
most likely correspond to the symmetry of (the critical 
points of) projection of the actual currents onto a surface 
determined by the collision axis, the collision plane, etc. 
These patterns well may have very little to do with the 
internal dynamics of the stable waveform. 



When mr 1, we generally have sinT — > and, ac- 
cording to Eq. (|2.10p , the pseudoscalar density nearly van- 
ishes and 5kK. Magnetic polarization of the Dirac field 
is greater than the electric one, \L\ > \K\. Vice versa, 
at the shortest possible distances, when T — > n/2, the 
Lorentz boosts play a major role. Accordingly, the pseu- 
doscalar density V ~ TZ is large and electric polarization 
is dominant, K 2 — L 2 ps 7Z 2 3> 1. These two sectors are 
separated by the two-dimensional surface T[i/j] = 7r/4. 
When T[tp] < 7r/4 (large distances) an appropriate choice 



for and will be vectors Hi and Hi of Eqs. (|2.7p . 

At the interval between A and V2X, where the mass of 

* 

the particle is being formed, these will be Ei and Ei, de- 
fined by a system of boosts, which are large along with 
the pseudoscalar density. 

In fact, this domain is also responsible for the magnetic 
moment of the Dirac particle. Indeed, using the classical 
formula for the magnetic moment, [i = (e/2c)[r x v], with 
the limit values, r = A = h/mc and v = c, which follow 
from the Eqs. (|4.29[) . we obtain 



' 2c mc 2mc 



(4.30) 



This is a well-known result corresponding to the gyro- 
magnetic ratio g = 2. 

The impossibility to introduce a normal orthogonal co- 
ordinate system in the presence of the axial potential N 
in the equations of motion is explicitly illustrated in Ap- 
pendix [Cj An attempt to separate the angular variables 
in the Dirac equation in the presence of only the radial 
component K r (r) is made. In this case, the radial coor- 
dinate, r, is a well defined normal coordinate. It appears 
that even in such a simplest case there are no opera- 
tors with eigenvalues of angular momentum that com- 
mute with the Hamiltonian. At the same time, angular 
variables can be explicitly separated in the equations of 
motion. The only possible explanation of this fact is that 
the formally introduced angles do not represent arcs of 
the usual spatial angular coordinates. There is no reason 
to require that solutions of the Dirac equations must be 
single-valued along these arcs. However, equations (|C.8[) 
for angular functions are clearly the equations for spher- 
ical harmonics. These harmonics can be interpreted only 
as elements of an internal symmetry in the space of po- 
larizations of the Dirac field. In order to find out what 
can be the non-geometric integrals of motion the angular 
and radial functions must be studied together. 



V. NONLINEAR DIRAC EQUATION. 

In this section, following the programme outlined in 
Scc lII B[ we will incorporate the nonlinear effects, which 
so far were found as constraints, into the Dirac equation. 
Following Fock [3], let us rewrite the operator a a tt a as 



a b n b 



(l/2)u> aca a c - (i/4:)e acb dOJ acb p 3 a d . (5.1) 
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Then, the Dirac equation reads as 



_d_ 

ds b 



+ ieA b + igp 3 ^b ~ \u aba 



-w b p 3 ip = -imp!ip, 
(5.2) 

where Hf, and w b = —(l/2)e ac db^acd are the sets of in- 
variants; the latter differ from zero whenever spacetime 
does not admit a coordinate net of all normal congru- 
ences. In general, the invariants w b do not vanish and 
they are complementary to the invariants H b given by 
Eqs. (|4.2ip and (14.22[) . [Since congruence is normal, 

we have w\ = -w 2 30, w 2 = —U310, w 3 = -uj 12 q.} 

According to Eqs. (|4.7|) and (|4.8p . each sum 
^2 a T}( a )Uj aba includes either the terms woio = dhiTZ/dsi 



(i = 1,2,3) or Ylfi^iM = dlnlZ/dso, all of which are 
bilinear forms of 4> + and ip and, geometrically, are the 
geodesic curvatures. At first glance, the presence of 
these curvatures makes the Dirac equation extremely 
non-linear. This genuine nonlinearity, however, can be 
effectively alleviated after the following "normalization" 
of the Dirac wave function. If we observe that 



dip 1 Shift 
ds a 2 ds a 



d_ 

ds a 



t 



(5.3) 



and assume Eq. (|4.17p be true, then we arrive at a 
much simpler equation for the normalized function £ = 
V'/v / ft=(3oo[V']) 1/ V: 



ds 



eAn 



2^3 [wo 



d a T] 



3 

£■ 

i=l 



d k 1 
i-gj- + i-^ - eAi + -pziwi + diT] 



mpi 



(5.4) 



where h 



E 



The nonlinear Dirac equation (|5.2|) now looks like a linear equation for the normalized 



function £. Once m -1 is accepted as a measure of length, this equation is dimensionless and does not change under 
a similarity transformation. At the first glance, the term p 3 diT in it is always nonlinear; however, the constraint 
k 3 = 2m sin T = 2m(£ + P2C) (cf. Eq. (|4.16[l ) eliminates the non-linearity whenever the curvature k 3 can be determined 
as a function of coordinates from geometric considerations. In the one body problem this curvature always has the 
meaning of the inverse radius of the enveloping convex surface. At least at large distances (at the scale of 1/m) 
we have T tx 1/mr and d r T oc — 1/mr 2 , which brings in a singular potential oc 1/r 2 into the Dirac equation and 
a Newton's force into Eqs. (|3.13p - (|3.16p . Such a construct is equivalent to the Newton's approximation in general 
relativity and it may serve as the first step of an iterative procedure for the Dirac equation. 

In general, there is no direct connection between the radius of curvature and the distance to any distinct point 
inside a localized object; the gradients diT can be arbitrary large. Most likely, the solutions of Eq. (|5.4p have multiple 
caustics where the invariant density 1Z is large and V dominates to the extent that 1Z ~ V. So far, we did not find in 
mathematical literature regular methods to study equations like (|5.4p . The methods of contact geometry [l4j seem 
to be most relevant. 

The most important source of the nonlinearity of the Dirac equation resides in that fact that evolution, in terms 
of proper time sq, has a different rate at different points of the localized object; the rate of evolution, d/dso, along 
this time cannot yield anything like the energy of this object as a whole. Fortunately, we have proved that a stable 
object does have a well defined hypersurface of a constant world time t. Therefore, a meaningful evolution scale for 
a stable object as a whole is associated with the physical time t. According to Eq. (|4.1ip . we have alt = IZdsQ. Hence, 



ill 



d_ 

di 



1 



eA ~ 7iP3[WQ + d T] 



3 

E' 



,d_ 

dsi 



1 



eAi + 7;PA w i + diT] 



mpi 



4> 



K 



(5.5) 



and now this equation has a scale fixing factor TZ in front 
of the d/dt. This effect is a major one because it corre- 
sponds to a clearly understood physical effect, refraction 
of the Dirac waves. It is of the same physical origin as 
self-focusing in nonlinear optics and acoustics or deflec- 
tion of light in the gravitational field of a star. In fact, TZ 
plays a role of "refractive index" depending on the ampli- 
tude. Since the phase velocity decreases with increasing 
amplitude, the field tends to auto-localize. This mech- 
anism of concentration is the most distinctive property 
of gravity which may signal its role in matter formation 



from fields at all scales. The spatially uniform solutions 
of the Dirac equation just cannot be stable. 



VI. LOCALIZED DIRAC MATTER, 
ELECTRIC CHARGE AND RADIATION. 

Our major perception regarding the vacuum is the ab- 
sence of matter. Since matter inevitably is localized, this 
means that in the vacuum, TZ is constant and the space- 
time metric of the Lorentz vacuum has goo = c 2 ~1. At 
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the same time, we have Eq. (|4.1ip , 

dt = 1Zds = dso/V5oo • P~TTj) 

Therefore, the empirical goo = 1 of an empty space cor- 
responds to 1Z = 1. This state cannot be stable. Due to 
a very special nonlinearity of Eq. (|5.5p Dirac waves tend 
to refract towards domains where TZ — 1 > amplify- 
ing TZ there until some saturation level (or caustic) is 
reached and an external boundary is formed. The oppo- 
site trend must be observed in domains where 1Z — 1 < 0; 
the Dirac waves tend to escape them. This conjecture 
can be phrased more precisely as: 

Identification of the sign of (1Z — 1) with the sign of elec- 
tric charge leads to a dynamic picture of an empirically 
known charge-asymmetric world in which stable positively 
charged on-mass-shell Dirac objects are highly localized 
(and presumably heavy) while negatively charged objects 
tend to be poorly localized (and presumably light). 
The best prospect of this idea is that these objects are 
the protons (or nuclei) and the electrons of the real world. 
When electric forces come into play, the electrons become 
somewhat localized around heavy objects, thus forming 
electrically neutral matter. For atomic electrons, which 
are smeared over distances much exceeding the Compton 
length and held weakly localized near nuclei by electric 
forces, the effect of 1Z — 1 <C 1 on the metric must be 
negligible. However, the view of a vacuum as a classical 
Dirac field with the standard level 1Z = 1 and propa- 
gating in it waveforms, instead of an ensemble of quan- 
tum oscillators with an unbound spectrum that interact 
with plane waves, may alleviate the problem of ultravi- 
olet divergences of the standard perturbation series and 
its renormalization. 

Our major conclusion about the nature of localization 
is drawn from the analysis of the on-mass-shell wave- 
forms. However, it relies on basic properties of the wave 
propagation so that it seems reasonable to apply them 
to at least the long-lived waveforms. From this perspec- 
tive, any positively charged particle should have a slightly 
longer lifetime and be more localized than its negative 
partner. While the proton is small and stable, the anti- 
proton should not have an as well defined outward bound- 
ary as the proton has. The lifetime of the anti-hydrogen 
may not be long even when it is completely isolated from 
normal matter. We refrain here from speculating about 
the possible mode of its decay. 

The next important question is the interaction of the 
localized Dirac waveforms with the electromagnetic field. 
Empirically, the field A c in the connection (|A.3|) is the 
electromagnetic field, which is responsible for the Lorentz 
force. The system of invariants H a in r a is determined 
either by geometric properties of congruences like (|4. 16|) 
or by nonlinear constraints (|4.2ip and (|4.22p . As long 
as the field H a is a gradient, the r.h.s. of Eq. (|3.9p is a 
system of invariants, which include, except for the geo- 
metric terms, the term ej a F a b. In the coordinate rep- 
resentation (|3.12|) . this part is translated into ej^F^, 
and this is the only term on the r.h.s. of the real part 



of Eq. (|3.9p . This term is known as the Lorentz force of 
the field F^ v — V^,, — V„A M acting on a charged par- 
ticle. The vector field = e^, A a originates from the 
connection T a and thus is an external field. 

For the real part of the energy-momentum tensor, 

Tab = {i/2)[^p + a a T)i > ip — tp + F) b a a, 0], and in the artificial 
normal coordinates, discussed in Sec lIIH we had equation 

±-y-g T %\=eJ-gfF, v . (6.1) 

As it should be, the kinematic acceleration = 
e (o)^ ?l,e (o) °f a criar g e d particle does not include a grav- 
itational part (see [lH, §63). 

The field F^ v is a tensor and it satisfies the identity 
(the first couple of Maxwell equations), 

F/xav = V M F W + V„F VI1 + V„iV = 0. (6.2) 

The divergence of this tensor, V^F^^ = 0, can be cast 
in the following form, 

— HF 4- R K F — R K F — R F Ka 

= (V^[V CT F TO ] - V^V'i^]) , (6.3) 

where V CT = g aX V \. The l.h.s. of this equation is the 
wave operator for the field F^ v . The two terms in the 
r.h.s. can be transformed further after we postulate the 
second couple of Maxwell equations, 

V M F^ = er , V^J" = 0, (6.4) 

with the Dirac's vector current eJ^ in the r.h.s. This 
amounts to the second (in fact, independent) definition 
of electric charge as the divergence of the electric field, 
and a few reservations must be made. First, without a 
good reason, the same coupling constant, as in the con- 
nection r a , is postulated. Only gauge invariance as an 
independent principle can provide for an unquestionable 
equality of these constants. Second, the Dirac field in 
Eq. (|6.3p is assumed to be a stable configuration and the 
field H M is considered a gradient. Otherwise, the con- 
servation of the vector current and the non-conservation 
of the axial current will conflict with Maxwell equations. 
Third, the interaction between the electromagnetic field 
and the spacelike axial current or pseudoscalar density, 
which are present in Eq. (|B.10p . and affect the balance of 
momenta inside the Dirac object, are disregarded. Only 
those interactions that are responsible for the change of 
timelike components of the momentum of a stable par- 
ticle are allowed to be sources of electromagnetic field. 
Then Eq.jlJl) becomes, 

— r\F 4- R K F — R K F — R F Ka — eO 

Q„u = (V M J v - V„ J M ), (6.5) 

where is a convenient intermediate notation. 

The crucial question is, if J ^ in this equation is (or 
can be) the current j^ of Eq. (|6.ip . which was derived 
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as a consequence of the Dirac equation with the poten- 
tial A c in the connection T c . Evidently, the answer is 
no because then, according to Eq. (|4.ip . we must have 
Qp,u = 0. Therefore, an on-mass-shell object cannot be 
a source of an electromagnetic field that may result in 
the Lorentz force of self-interaction. The definition (|6.4j) 
must be complemented by Eq. (|6.1[) . which then deter- 
mines the measured acceleration of another charge that 
senses the field of the first one. The potential A c in the 
Dirac equation that describes a localized object must be 
"external" ; its source can be only the current of another 
object. The problems of mass and charge (including 
the problem of electromagnetic mass) are not a single- 
body problem. One further implication of this observa- 
tion is that if one can simultaneously identify two local- 
ized objects, then the Dirac fields of these objects can- 
not overlap in space-time . Since, for a stable object, 
one cannot set ej v — V^F''" in the expression for the 
Lorentz force, it is also impossible to express this force 
as the divergence of the energy-momentum tensor, i.e., as 
ej^F^ = V^Fp, = -V^e^F) and claim that 6° 
is the energy density of the electromagnetic field. This is 
not surprising since one cannot convert this energy into 
any other form without a second object. 

For an isolated stable charged object the condition 
that it cannot interact with its own electric field means 
that the only "potential" in the wave equation (|5.2[) is 
N oc 1/r 2 . The wave equation with such a steep potential 
may have a strongly localized solution regardless of the 
sign of this potential. The difference from the commonly 
studied cases is that now the signs of this potential for the 
left- and right-handed components are opposite and that, 
for a stable wave form, the region mr < 1 is cut off by 
Eq. (|4.16[) as unphysical. Therefore, a precursor of a lo- 
calized state is present in the Dirac equation even before 
the universal nonlinear mechanism of the time slowdown 
takes over. 

The field which is measured via the Lorentz force (|6.1|) 
always is a "field in vacuum" . The wave equation (|6.5|) 
for the F^u from Eq. (|6.ip is a homogeneous equation, 
which depends on the Dirac field of (|6.ip only paramet- 
rically, via derivatives of (^(VO i n the Ricmann tensor. 
The electromagnetic sector of the theory turns out to be 
entirely in the form required by Riemannian geometry. 
This sector is responsible for the propagation of signals 
that are used to synchronize macroscopic clocks (which is 
unambiguous only in special relativity). 4 A stable wave- 
form of the Dirac field with — neither interacts 



3 From perspective of the second quantization, when eigenfunc- 
tions of the Dirac equation are associated with different states, 
this means that the Fock operators of these two states must anti- 
commute! 

4 It is important to emphasize that the l.h.s. of Eq. ll6.5t is given 
in terms of measurable electric and magnetic fields; therefore we 
indeed are dealing with a signal that may have leading and rear 
fronts. Since the world time T 114.101 1 is a harmonic function it 
can be discontinuous along characteristics. The time of emission 



with its own electromagnetic field nor can it emit an 
electromagnetic field, as a signal, by itself. This is yet 
further evidence that the object is in a stationary state. 
This property is in line with the well known fact that 
equation of the Coulomb's law is a constraint and not an 
equation of motion. The longitudinal part of the elec- 
tromagnetic field does not propagate; the Coulomb field 
is simultaneous with its source. The field of radiation 
emerges only when this simultaneity is lost. This is yet 
another view of the realm of the well-known phenomena 
of transient processes where the proper field of a particle 
is truncated [ill]. 

What if it occurred possible to trace an observed ra- 
diation field back to the current in the interior of the 
localized object (so that = j^) and, e.g. by a pre- 
cise analysis of radiation, to learn that 7^ there? 
Then Eq. (|4.1[) must be replaced by the second equation 
of (|6.5[) . Proceeding as previously, we get 

V^eW - V„eW + e£H In ft - efft, In TZ 

= -(l/K)Q^. (6.6) 

Contracting this equation with spacelike e (i) e Q,) ih3 = 
1, 2, 3) and recalling Eqs. (|2~T5|) we find that 

woy - ojQji = —(l/TZ)Qij . (6.7) 

If Qij starting from some time moment to, then at 
t > to the congruence ety, of lines of the vector current 
cannot be a normal congruence. The family of space- 
like surfaces t = const, orthogonal to the vector current, 
vanishes. This means that Eq. (|4.4[) cannot be obtained 
and Dirac field cannot form a stable object 5 . The elec- 
tromagnetic fields produced by such an object are not 
just longitudinal (Coulomb) fields and the object must 
start to radiate solely because the electromagnetic field 
around it is not simultaneous with its source. Since the 
Dirac equation is of the hyperbolic type, the changes of 
the Dirac field must propagate also, having a light cone 
as a leading wave front. 

Contracting Eq. (|6.6[) with e^e^ we obtain another 
equation, 

Wioo = -{dkiK/dsi) - (l/n)Q i0 , (6.8) 



of a photon, in principle, is not defined. A photon does not 
constitute a signal. 

This is yet another way to view two seemingly different phenom- 
ena, the Meissner effect and precession of the spin in magnetic 
field. In order to be in a stable quantum state, the supercon- 
ductor expels magnetic field from its interior or confines it into 
vortices, thus defining a common time across its whole volume. 
In the same way, electron with the magnetic moment 1 14.301 , be- 
ing placed in magnetic field, moves in precession with the Larmor 
frequency. Therefore, in rotating frame the magnetic field van- 
ishes and the electron still can have the same world time across 
its volume (staying in a certain quantum state). 
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that accounts for the effect of the electric field, which 
adds a boost in the direction of the congruence . In- 
teraction with the electric field alone (which can be the 
case only when this field is longitudinal) does not destroy 
the hypersurfaces of constant time of a localized object, 
which allows it to stay intact. The most important effect 
of acceleration in an electric field is altering the shape of 
a charged object which leads to an increase of its internal 
energy and of the local charge density. 

Referring to the above qualitative analysis we may go 
further and discuss a qualitative picture of some transient 
processes. If the Dirac waveform is not stable (as is in 
the case of then the development of instability (and 
the lifetime) must still be stretched, due to the nonlinear 
effect of the time slowdown. When the limit of stability 
(bifurcation) at r ~ Aju = 1.86 • 10~ 13 cm is reached, then 
the previous dynamic regime suddenly breaks up and the 
field begins to evolve towards a new configuration of a 
smaller mass m e and a larger A e = 3.86 • 10 _11 cm. The 
Dirac field of the which originally occupied the in- 
terval A M < r < A e , must be radiated. Since the Dirac 
equation is hyperbolic, the sharp front of the Dirac field, 
as any precursor, must propagate along a characteristic 
(at the speed of light), in the outward direction. This can 
only be the right component of the Dirac spinor (which 
inherits its lightlike current from the caustic), the v^. 
The final state of a e + emerges not earlier than a new 
caustic at r — A e is formed. This requires yet another 
transient process of the violent collapse of the Dirac mat- 
ter onto a new caustic and radiation of a similar pre- 
cursor, but with the opposite phase velocity, v e . The 
leptonic number is preserved dynamically in both pro- 
cesses. Remarkably, it is exactly the existence of a well- 
defined (by the spacelike axial vector) outward direction 
that eliminates the illusion of the reflection symmetry 
of a plane wave and thus predetermines a unique po- 
larization of the spinor precursors. As it was noted by 
Wigner [UJ , it is only a theoretical idea of mirror symme- 
try (expressed in terms of polar vectors) that hints of the 
possible existence of the second polarization for lightlikc 
spinor waveforms. 



VII. CONCLUSION. 

The nonlinear Dirac equation, with its capricious inter- 
play of the many polarization degrees of freedom, poses 
a tough mathematical challenge for theory. Its explicit 
solutions may well yield various "magic numbers" that 
are currently known only from experiment. Even before 
regular mathematical methods are developed, one may 
rely on various qualitative consequences of the finite size 
of the Dirac waveforms to re-analyze existing data. 

1. The conjectured connection between the mech- 
anism of self-localization and the sign of the electric 
charge of the Dirac wave form also assumes that posi- 
tively charged particles, which are not perfectly stable, 
must have a somewhat longer lifetime than their nega- 



tively charged anti-particles. The ratio \ — (t + — r_)/r a „ 
was measured for the most long-lived species as a test 
of CPT-invariance. According to the Particle Data 
Group, the difference in lifetime is indeed always positive, 
X{K±) = (0.11 ± 0.09) • 10~ 2 , x(^) = (5-5 ± 7.1) • 10~ 4 , 
x(/-* ) = (2± 11) • 10 -5 , being the largest for the heaviest 
specie. 

2. One of the predicted manifestations of charge asym- 
metry is the existence of the particle's external size. The 
internal radius is universally limited from the below by 
the Compton length A = h/mc. Due to the time slow- 
down in domains of large Dirac density, the positively 
charged species must have smaller external size than their 
negatively charged partners. Possibly, e + has reasonably 
well defined external boundary, which then may explain 
its relatively long lifetime in the environment of normal 
matter. There may well exist observed differences in the 
dynamics of electrons and positrons (or p and p, e.g., in 
accelerators) that are currently attributed to technical 
issues. 

3. There are certain coincidences of numbers that may 
prompt another look at well known phenomena. We 
know that, A e = 386/m, A M = 1.86/m. The radius of the 
proton, as estimated via its electromagnetic form-factors, 
is r p = lfm and it grows for nuclei as A 1 / 3 . This may 
well tell us something about the high rate of \i~ capture 
by light nuclei versus the low rate of the inverse /?-decay 
by even heavy nuclei. The correlation between the cap- 
ture rate and size of a nucleus may be a useful test. 

4. The failure to keep anti-hydrogen molecule in the 
cold atom trap for an indefinitely long time may establish 
the limits of stability of the anti-proton in anti-matter 
surroundings. One may think of the capture of e + by a 
poorly localized p as a possible mechanisms of instability. 

5. The Lorentz contraction of the accelerated wave- 
forms is a dynamic effect, which leads to the accumu- 
lation of energy in an extremely small volume and its 
release in the course of a collision. The 5-8 Gev electrons 
and positrons are compressed to a size about 10 /m 
having a density higher than a proton. Upon colliding, 
they stop and create a sharp peak of invariant density, 
which is very far from a stable configuration and rapidly 
decays. The B ± lifetime is reasonably long, 1.7 • 10~ 12 s, 
and its size is about 10 _1 /m. The time slowdown at 
1Z ^> 1 may result in two effects: (i) an abnormally long 
lifetime of the resonance and (ii) an exotic trend to fur- 
ther decay into compact heavy objects rather than to 
decay into lighter objects according to the usual specta- 
tor model. The mode B+ -> A+A(3872) -> J/ipir+ir- 
seems to be a candidate for this kind of the process be- 
cause the width of A (3872) is very small. 



APPENDIX A: PARALLEL TRANSPORT OF 
DIRAC FIELD. 

In order to derive a measure for comparison of the 
fields ip{x\) and ipix?) at two close points let us require, 
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following Fock 0], that the components j a = ip + a a ip 
are the invariants of the vector j^{x) and the con- 
gruences e? a \(x) at the point where vector is defined, 

j a = e^j 11 . For now, we assume that a set of four or- 
thogonal congruences is fixed in advance and that Dirac 
matrices a a are either invariants or covariantly constant 
objects. When the invariant j a is parallel-transported 
by dsb along an arc of congruence (6), then, solely be- 
cause the local pyramid is being rotated, it must change 
by 5j a = oJ a cbj c ds b = uj acb ip + a c i/jds b . The invari- 
ants ujabc are defined by Eq. ()2.14p . Let matrix T a (the 
connection) define the change of the Dirac field compo- 
nents in the course of the same infinitesimal displace- 
ment, Sip — T a ?pds a , 5%jj + = ijj + T+ds a . Let differential 
of the product ip + a°"ip obey Leibnitz rule. This gives yet 
another expression for 8j a , 

S Ja = ^ + {T+a a + aa T b )4>ds b . (A.l) 

The two forms of Sj a must be identical. Hence, the equa- 
tion that defines T a is 

r+a Q + a a T b = u) acb a c , (A. 2) 

and it has the most general solution, 

T b (x) = ieA b (x) + igp 3 H b (x) 
1 i 

--^OJ0kb(x)p3^k ~ -^Okim^irnb{x)(Jk , (A. 3) 

where the last two terms can be compacted as, fl b = 
(\/<i)uj cdb pia c pia d = (l/4)uj cdb j c j d . These two terms 
correspond to an infinitesimal boost of (|2.4[) along the 
spatial fc-axis with parameter uj k b ds b and an infinitesi- 
mal rotation of (|2.3[) in the (im)-plane with parameter 
Wi m bds b , respectively. This analogy, however, is limited. 
While Eqs. (|A.ip - (|A.3[) do imply some measure for the 
length of an arc (and of an angle as the ratio of the two 
lengths), Eq. (|2.2p does not. The first two terms are due 
to an intrinsic indeterminacy that arises when one has 
to compare Dirac fields at two different points relying 
only on the properties of the vector forms (ip) . The 
first term is readily associated with the electromagnetic 
potential. The second one would not appear at all if, 
following Fock 0], we required that S(tp + pi^p) = and 
5(ijj + P21P) = 0- This decision was motivated by that kind 
of invariance of the Dirac equation in Minkowski space 
(local Lorentz invariance), which is not inherited by the 
Dirac field in Riemannian geometry. The position of 
in connection (|A.3[) may lead to the impression that it 
can well be a "next field" , which interacts with the ax- 
ial current J b of the Dirac field and is governed by an 
independent equation of motion. At least for the on- 
mass-shell Dirac field, this is not true. 

The connection (|A.3|) commutes with the matrix p% 
so that Eq. (|A.2[) remains the same when a a — > pza a - It 
neither commutes nor anti-commutes with p\ and p2, viz. 

r b>i + Pi^b = 2gp 2 K b , 

r fc + p 2 + P2 r b = -2g Pl X b . (A.4) 



From now on, we postulate that invariant derivative of 
the Dirac field is D a ip = (d a — T a )-0 where d a — e^<9 M is 
the derivative in the direction of a curve of congruence 
(a). Assuming the Leibnitz rule for D a and considering 
all Dirac matrices as constants we readily reproduce the 
reference point of Eqs. (|A.ip and (|A.2|) as 

D b j a = d b j a - UJacbjc = ^bja ■ (A. 5) 

The result (|A.5|) for D b j a is a warrant that after pro- 
jecting the r.h.s. into coordinate space we must re- 
cover the covariant derivative V^jt, . Indeed, e a e y V b j a = 
dftjv — r^Jer = V ' nj v , and we shall consider this as a 
proof that j a is an invariant of the vector and con- 
gruence e^) 6 - (There is no one-to-one correspondence 
between the two terms of V b j a and V^j,,!) 

In exactly the same way we may verify that the invari- 
ants D b J a of the axial current are of the form T) b J a = 
^ b J a and conclude that e^e^DhJa = V ' pj v . Vice versa, 
the quantities D b J a = e^'e^V y,J v are invariants of a ten- 
sor and a system of congruences. It is straightforward to 
verify (computing all derivatives as functions of ip) that 
equations like (|A.5|) hold not only for and but 
for e^[V>] = j^/TZ and e^ 3 J^] = J* 1 /1Z. Recalling dis- 
cussion of Sec. II, we may view this as complementary to 
(I2.13[) . i.e., proof that the unit vectors e^L and are 
vectors of Riemannian geometry. For the same reason, 
the projectors [5£ - fj v jV?\ and [6% + J^J„/K 2 } are 
the tensors. 

Using the same technique of differentiating and by 
virtue of Eqs. (|A.4[) we obtain 

D a S = d a S - 2gVK, D a V = 8 a V + 2gSK- (A.6) 

As one can see, that there is no immediate correspon- 
dence between the algebraic and differential properties 
of the scalars. However, the quantities from the first 
line of (|2.8[) (like TV) are differentiated as true scalars. 
The same behavior is observed for the components of 
the skew-symmetric tensor M ab . Instead of the antici- 
pated e c x e^ 1 e b J D c M a b = Va-M^j, we encounter one more 
disagreement with the differential criterion (|A.5|) . 

DcMab = ^cMab - 2gK c Mab, (A.7) 
D c Mab = VcMab + 2gtt c Mab- 

We leave open the question of if and when and of 
Eqs. (|2.7p are the vectors with the same degree of confi- 
dence as and e^ 3 %. In most cases, and e^ 2 % corre- 
spond to angular coordinates, which are physically uncer- 
tain without external fields (other objects nearby). The 
additional terms in Eqs. (|A.7p indicate that these congru- 
ences can be normal only in very rare cases of solutions 



Schouten (0, Ch.II, §9; Ch.III, §9) considers equations like (lA75l 
as a condition that fixes the components of a vector with respect 
to anholonomic coordinate system. 
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with a dynamically generated symmetry. The expected 
one-to-one match with geometry is spoiled by the extra 
(with respect to generator of Lorentz transformations in 
the connection r a ) matrices p\ and p 2 responsible for 
the "mixing" between right and left components. From 
this perspective, the Sakharov's idea Q, regarding the 
topological nature of elementary charges, seems be closer 
to reality that it initially appeared. A further analysis 
of this issue seems to be impossible without an explicit 
solution of the equations of motion. 

Recalling Eqs. (|2.9[) and using (|A.6|) we can compare 
the covariant derivative D a V computed in two ways, as 
D a V = d a {7l ■ sinT) + 2<?7£cosYN a or, alternatively, as 
D a V = d a lZ ■ sin T + 1Z cos TD a T. The result reads as 

D a T[ip] = d a T + 2gK, (A.8) 

which is invariant under the simultaneous transforma- 
tions, T -► T + Y{x) and 2<?tt a -> 2 S H a - d a Y(x). A 
supposed freedom of such (chiral) transformations is not 
permissible, since these transformations change the ob- 
servables in the r.h.s. of Eqs. (|3.13|) and (|B.4[) without 
altering the l.h.s. 

The commutator [D a , Z3 b ] still contains derivatives. In- 
deed, 

[D a , D b ]ijj = {d a d b - <9 b <9 Q )-0 

-[9 a r fc -9 b r a ~r a r h + r fc r a ]^. (A.9) 

However, for practical purposes it is important that 
[D a ,D b ]ijj can be split into two parts, with and with- 
out derivatives. Since ip is a coordinate scalar and, in 



general, derivatives along arcs do not commute, we have 
{d a d b - d b d a )ip = {uj cab - co cba )d c ip. Now we can re- 
assemble [D a ,D b ]ip as follows, 

a , l3 b ]4> = {Uj cab - LU cba )F) c Tp - D ab , (A.10) 

o ab = [d a v b - d b v a - r a r b + T b r a - c cab r c }^, 

where the term C cab T c = {uj cab — uj cba )T c is added and 
subtracted to replace d c tp by the covariant derivative 
D c tfj. The matrix operator B ab = — e^e^D^, D v ] in the 
r.h.s. does not contain derivatives and can be explicitly 
found, 

Dab = -^RabcdPiof pia d + ieF ab + igp 3 U ab , (A. 11) 

F a b = d a A b — d b A a — (uj cab — Lo cba )A c , 
Uab = <9 a K b - <9 b K a - (u cab - w cba )H c , 

where invariants of the Ricmann curvature tensor are de- 
fined by Eq. (|2.22p and F ab and U a b are invariants of the 
electromagnetic tensor F^ — V ^A^ — V ' V A^ and the field 
tensor — V^Kj, — Vj,K M , respectively. 

Using the cyclic symmetry of the Riemannian curva- 
ture tensor it is straightforward to show 0] that 

a a D ab = ^a a R ab + iea a F ab + igp 3 a a U ab . (A.12) 

When field H Q in the connection T a is a gradient, the ten- 
sor U^v vanishes identically, which is assumed throughout 
this paper except for Eq. (|B . 10|) , where a possible connec- 
tion with electro-weak theory is mentioned. 



I 

APPENDIX B: STRESS TENSOR AND PION FIELD. 

1. Internal flux of mass and stress in the Dirac field. 

In this section we study the stress tensor P b = iip + p 3 a a D b ip, mostly following the same logic as for the energy 
momentum tensor T b a = iip + a a D b ip in Sec lIIIB . starting from its covariant derivative. We find that 

D c [^ + p 3 a a ^ b ^] = d c [iP + a a D* b i>] - cj adc ^+ p 3 a d ^. (B.l) 

Once again, the last term of Eq. (|3.5p is missing, and thus we have no confidence that the covariant derivative is a 
tensor. This time, let us begin by contracting indices a and b in Eq. (|B.l( ). 

D c [^ + p 3 a a D* a ^] = d c [i/j + p 3 a a D* a i>} - Lu abc ij + p 3 a b l5 a ip. (B.2) 

By virtue of the Dirac equations, the first term in the r.h.s. of (1B.2|I becomes d c [mtp + p 2 4'}- Alternatively, we can 
immediately use the equations of motion in the l.h.s. and only then differentiate (matrices p 3 and a a commute), 

D c [^ + p 3 a a D* a i>} = mD^+p^j] = md c [^ + p 2 ^] + m ■ 2 5 5H C . (B.3) 
Comparing the last two equations we finally get the equation, 

uJacb ■ P ca = -2igmS$ b , (B.4) 

which is complementary to Eq. (|3.13[) . The imaginary part in the l.h.s. is due to (l/2)[P ca — P+J = (i/2)D c J a . Since 
the axial current is a vector, we can rewrite the last equation as 



(l/2)w Qcb V c J Q = -2gmSX b , 



(B.5) 
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which is complementary (dual) to Eq. (|4.20p . The skew-symmetric Hermitian part, (P ca + P c +) — (P ac + P^ c ), must 
vanish since the r.h.s. of Eq. (|B.4[) is an imaginary quantity. This yields the equation, which duplicates Eq. (|4.1[) . 



i[ip + p 3 a a l3 c ip - ip +< D^ a a p 3 ->p - tp + p 3 a c l5 ' a ip + ^ + ^D^a c p 3 ijj] = e acut D u j t = 0. 

and thus indicates that we still are dealing with a stable waveform. 

Contracting (in Eq. (|B.2p ) indices a and c we arrive at the expression, which is similar to Eq.([3 



D a [^+p 3 a a ^ b 4,} = d a [iP+p 3 a a ^ b i/j}+Lu acc yj+p 3 a a Tl b 4, 
Let us first rewrite the l.h.s. of Eq. (1B.7|) as 



1 d 



-g dx v 



(B.6) 



(B.7) 



D a [^+p 3 a a ^ b ^} = i/j+p 3 a a a ^ b - ~d b ~d a ]il) + 4>+*D + a p 3 a a ^ b ^ + D b {^+ p 3 a a ^ - ^ +< D + b p 3 a a ^ : (B.8) 

Because of an obvious change of signs (caused by an extra p 3 ), the last three terms of (|B.8[) do not cancel. Instead of 
(|3.9P we have 

D a P b a = i^+p 3 a a a ^ b - T5 b T5 a ]^ + imD b T + tm[^ + p 2 '^ b ^ - P2# (B.9) 

By splitting the commutator according to (|A.10[) . we can assemble the covariant derivative in the l.h.s. as 

= V M Re(P£) + l -W^ u J^ = V M Re(P£) + % -J»R^ + imd„V, 

leaving on the r.h.s a remainder, e^[w ca6 P c a + eJ a F ab + gj a U ab + {i/2)J a R ab \. By virtue of Eqs. (IB~4|) and (|X6|) the 
imaginary terms on both sides exactly cancel each other and the remaining real part reads as 



V M Re(P^) = ejW^ + gfU„ v + im[^p 2 ~6 v ^ - V; + S+p 2 # 



(B.10) 



The flux of momentum in the spacelike direction is de- 
termined by the Lorentz force that acts on the electric 
axial current ej^, but also by the Lorentz force of the 
field that acts on the vector current gj^. This part 
is logical to write down in terms of the mixed fields that 
act on lightlike left and right Dirac currents. It describes 
the interplay between the right and left components of 
a Dirac particle and belongs to the area of electro-weak 
interactions. Being interested only in stable states, we do 
not consider it in any details here. The last term is due 
to convection transport of the pseudoscalar mass density 
rnP in an electromagnetic field. An importance of axial 
electric forces for the processes of pion electro-production 
was noticed by Nambu and Shrauner [l5|]. The convec- 
tion term can be cast as 

* cd 

m[^ + p 2 {id b ip) - (id b ip + )p 2 4> + 2eA b V + (l/2)uj cdb M }, 

which repeats the familiar pattern of Gordon's decom- 
position of the Dirac (vector) current with the replace- 
ment pi — > p2, e — s- m, and where the long derivative 
includes only the electromagnetic potential. The pseu- 
doscalar density V is one of many polarization degrees 
of freedom of the Dirac field and is not an independent 
field. In some approximation, the charged pseudoscalar 
flux inside Dirac waveforms (e.g., nuclei) can be viewed 
as the interaction of highly localized nucleons via soft 



pion exchange. Complementary positions of vector and 
axial currents in Eq. (|B.10[) prompt a parallel between the 
charge and chirality of the Dirac field. This parallel was 
a point of departure for the model of elementary particles 
developed by Nambu and Jona-Lasinio (l6j . 

Eq. (|B.10[) also describes a process within a compact 
object that take place during a period of its acceleration. 
In its course, the object changes configuration, undergoes 
Lorentz contraction and, in fact, becomes a different ob- 
ject, with a larger density and slower flowing time in its 
interior. Therefore, the annoying non-unitarity of the 
proper Lorentz transformations is a physical effect. 



2. Pion field and tv <-» 27 processes. 

The flux of charge, mass and momentum carried by 
the pseudoscalar density in the interaction between Dirac 
nucleons is commonly attributed to the pion field. Pions 
can also be detected as sufficiently long lived particles. 
Therefore, V should satisfy the Klein-Gordon equation, 
which we will derive below. 

The first step is to put the axial current in a form 
with separated convection and polarization currents, as 
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is done in Gordon's decomposition of the vector current, 

r = -^[-iP + a [a p 2 a b] L) b iP - ^ +< D^a [ V 2 a 6 V]- 

where o\ a ...a b ' i stands for a a ...a b — a b ...a a . (Now, the 
entire convection term is reduced to the derivative of VI) 
Computing the covariant derivative of both sides of the 
last equation and using (|3.4p we obtain 

D\V = 2^ +< D + p 2 ~S a 4> - 2m 2 V 
-Re[^+J a p 2 a b ^ (D a D b - D b D a )^]. (B.12) 

For the on-mass-shell Dirac field of the nucleon with 
a large mass to, we may take the Dirac field in semi- 
classical approximation, ip cx e lS ^ n , so that the first two 
terms in the r.h.s. constitute the classical Hamilton- 
Jacobi equation for the eikonal S. If this equation is sat- 
isfied with a sufficient accuracy (the waveform behaves as 
a classical particle and the resonance is sufficiently nar- 
row), then in the r.h.s. remains only the last term. By 
virtue of Eqs. (|A.10p and (|A. 1 1|) this term becomes noth- 
ing but eF ab M ab + (R S /2)V - (l/2)C cab ^+a a p 2 a b ^ c ^, 
where R s = R[4> + ,ip] is the scalar Riemannian curvature 
(with dimension to 2 ), which is a functional of the Dirac 
waveform. As a result, we arrive at the independent wave 
equation for the pseudoscalar density V (the pion field), 



D 2 a V - y-P R« -C cab Re[i; + a a p 2 a b ^ c ^} 



+eF ab M 



ah 



(B.13) 



a Heisenberg equation of motion with a variable mass 
defined by the negative scalar Riemannian curvature R s 
outside the stable nucleons 7 . Quite surprisingly, exactly 
this V enters the r.h.s. of Eq. (|3.18p that defines the force 
of gravity/inertia in the same approximation of a mate- 
rial point. 

The source in the r.h.s. is Hermitian. Its first term 
is "geometric" and accounts for the flux of momentum 
and twist of the tetrad basis which are needed in order 
that the pion could be born. The second term is more 
related to the pion decay and it can be rewritten in two 



7 The negative curvature is typical for the geometry of the expand- 
ing matter. The pions (and mostly pions) are abundantly cre- 
ated in the high-energy collisions where strong contraction of the 
colliding particles is translated into the rapidity plateau in the 
distribution of pions. At any given moment of time t, the proper 
time at a distance x from a generic point x = corresponds to 
the earlier proper time r and has a larger density 7Z(x, t). There- 
fore, as can be perceived from any point, the proper time flows 
more slowly with larger distance x from this point. As a result, 
Dirac waves tend to refract into distant spatial regions. The issue 
of stability of the expanding Dirac field and of its localization as 
pions and muons will be discussed elsewhere. 



ways. On the one hand it can be presented as 2etp + [piE+ 
p 2 B] ■ aip. When F ab is the field of a standing transverse 
electromagnetic wave it has a simple representation in 
terms of the spin interaction with two waves of circular 
polarizations, 



eF ab M ab = 4eJ2 



■■(C k e 



-ikx sy* Ahx 



x [e L ■ tp R dipL + e R ■ ip L aip R }, 



C* k e lkx ) 



(B.14) 



where eL,ii(k)-Lk are the vectors of the two circular po- 
larizations, ipL.R are the left and right components of the 
Dirac spinor field, lo^ = is the "photon's energy" and 
Cfe is the Fourier component of the initial or final (possi- 
bly, coherent) state of the Heisenberg field F ab . This form 
of the source of the pion field allows one to qualify ir° as a 
resonance in the system of the Dirac field and a standing 
electromagnetic wave formed by the two circular polar- 
izations, which causes the simultaneous flip of helicity of 
both components of the Dirac field. On the other hand 
the source in the wave equation (|B.13[) has the structure 
of the axial anomaly. This could be an exact correspon- 
dence if there was a simple proportionality between Ai ab 
and F ab . Then, eF ab M ab = CeF ab F ab , where the ex- 
plicit value of C must comply with the observed rate of 
the 7T° — > 27 decay. It is instructive that the wave equa- 
tion for the pseudoscalar meson field V (that yields the 
pole in the pion propagator) was derived exactly from 
the original equation Eq. (|3.4| . The term, which was ad 
hoc added to this equation by S. Adler pj| (in order to 
save the Ward identity for the axial vertex in triangle 
graph) has naturally appeared as the source in the wave 
equation. A uniform method to derive wave equations 
for meson fields is beyond the scope of this paper (2(| • 

In the first approximation, the value of mass of the 
Dirac field is not important. For this particular reso- 
nance, the mass term in the l.h.s. of Eq. (|B.13j) can be 
confidently identified with and measured as (2E 1 ) 2 . In 
fact, m 2 ~ — R s /2 is an independent of the Dirac mass 
to measure of the "metric elasticity" in the ground state 
of the Dirac field, when balance between left and right is 
probed by electromagnetic field, thus being a fundamen- 
tal constant. The geometry of currents inside 7r° as a 
finite-sized object is not clear so far. The dynamic quan- 
tity m v is meaningful only for 10~ 16 s of the resonance 
spike of the pseudoscalar density. The form of the source 
(sink) in the Eq. (|B.13p indicates that this spike is a para- 
metric resonance which can be created in a reaction like 
7P — > 7r°p — > 77P (or, more precisely, 7*7* — > 7r° — > 77). 
It decays due to axial polarization currents (eventually 
producing two photons of the opposite circular polariza- 
tion) . The nucleon is needed solely to trigger the process 
- to locally determine the difference between left and 
right or inward and outward. Another way to excite the 
7T° resonance is e + e _ -annihilation at high energy, etc. 
For the sake of completeness mention that the source in 
Eq. (|B.13p is complementary to the interaction of the ax- 
ial current with the electromagnetic field given by the 
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term eJ^F^ v in Eq. (|B.10p . which describes the balance 
of momentum when the waveform changes its shape. 

The processes with the charged pions can be consid- 
ered in a similar manner [201 ] . The dynamic footing for 
the phenomenological formalism of isospin symmetry, as 
it is presented in various texts (e.g. [ll|), can be naturally 
built from the basic properties of the Dirac field. Pions 
are identified as one specific (pseudoscalar) polarization 
degree of freedom of the Dirac waveforms. The totally 



dynamic origin of the pion mass term, which is deter- 
mined by the curvature R s [tp + , ip] of the matter- induced 
metric, explains the diversity of faces that pions may re- 
veal in different situations. This variety ranges from soft 
pion glue within nuclei (when DIS cannot resolve pion's 
structure functions) and up to free propagation of the 
massive pions (tracks) at distances that allow for the pion 
interferometry (which is sensitive to the microscopic dy- 
namics of the pions emission [22| ) . 



I 

APPENDIX C: ANGULAR VARIABLES IN DIRAC EQUATION. 

Let us examine the properties of the solutions of the Dirac equation (|5.3p . neglecting nonlinear terms, in the presence 
of the radial field 2gH r = — <9 r T and in a perfectly spherically symmetric geometry. Since we will focus on the nature 
of angular variables, an explicit dependence H r (r) is not essential. [One may think of Eq. (|4.18|) as an example.] The 
only non- vanishing components of the Ricci rotation coefficients are W212 = (l/ r ) c °t ^ an( i w i3i — ^232 = 1/r. Solely 
as a reference, assume that there is an external electromagnetic Coulomb field Ao(r). Then the Dirac equation is 

\id - eA + ga 3 K - ip3cr 3 (d r + -) - i^^-(d g + -cot0) - i Pza ' 1 d v - mpx\i/> = . (C.l) 
L r r 2 rsmO J 

In terms of a new unknown function, ip(r, 8, ip) = r\/sin Qip , this equation becomes 

lido ~ eA + gcr 3 H r + p 3 o- 3 (-id r ) + — (-ia-^de - i-^cL) - mpAip = . (C.2) 

L <r Gin H J 



Equation (|C.1[) is the Dirac equation in the tetrad ba- 
sis. In order to find its solution one has to separate the 
angular and radial variables. This is known to be a some- 
what tricky problem, even in the standard problem with 
a radial Coulomb field(when H a = ). The Hermitian 
operators in Eq. (|C.2|) are the tetrad components of the 



momenta j<3 = —id r , p\ = 



P2 



-i(r shri 



The operators p\ and P2 are clearly associated with the 
angular motion. If the coefficients in this equation where 
not matrices, it would have already been an equation 
with separated variables, which would match the per- 
fect spherical symmetry of the external fields A$ and 
N r . The problem is that the operators of the radial and 
angular momenta do not commute (they anti-commutc, 
[a3j?3, (aipi +a.2P2)]+ = 0). A regular way to avoid this 
obstacle is as follows pj], [H| : One attempts to construct 
a minimal set of operators that commute with the Hamil- 
tonian. For example, one can check that the commuta- 
tor [03^3, pi(ctipi + a2f>2)]- = and take the operator 
pi(a\pi + ot2P2) as a generator of the conserved quan- 
tum number. This trick works when H = and it is very 
instructive to see the details of its failure when H 7^ 0. 

The conventional operator of angular momentum is 
C = [fxp] +a/2. An additional operator C — a- C— 1/2 
commutes with the orbital momentum, [C, (r x p)\ = 0, 
and has the properties, £(£ — 1) = [r x p\ 2 and £ 2 = 
£ 2 + 1/4. Therefore, if n is an eigenvalue of operator £ 
we obviously have k(k — 1) = 1(1 + 1) and n 2 > 0. On 
the other hand, if Ca = Pi£, then (Ca) 2 = £ 2 and these 
operators have the same sets of eigenvalues. In the tetrad 



basis, these generators of the angular quantum numbers 



are 



£ 



Pi ( 



ia 2 d e 



l(J\ 



-id. 



^-a 3 .(C3) 



In terms of the auxiliary operator Ca (which has the same 
set of quantum numbers as the operator of the angular 
momentum but is a different operator) and the projection 
£3 of angular momentum, the Dirac equation becomes 



[ido - eA + .gcr 3 N r + pzo- 3 (-id r ) 
Ca 



-P203- 



mpi]ijj = 0. 



(C.4) 



When the operator Ca commutes with all terms of 
Hamiltonian (which is the case when H M = 0) we can re- 
quire the wave function be an eigenfunction of the Hamil- 
tonian and these two operators, 



CaiP — k1 Pj an d 



£# = K + \)^. (C.5) 



Since the Ca anti-commutes with (703 K r we have no ob- 
vious solution for the separation of variables. 

Because the presence of the component H r (r) at least 
apparently preserves the spherical symmetry, we can try 
to look for a general solution of the following form (£ and 
r] are the left and right components of the Dirac field in 
spinor representation, respectively), 



u L (r,t)y(9,cp) 
d L (r,t)Z(6,<p) 



u R (r,t)y(6,<pf 
d R (r,t)Z(9,<py 



(C.6) 
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As a first step, we may try to substitute the Dirac spinor 
(|C.6p into Eqs. (|C.5p . One can immediately see that the 
angular variables in (|C.5j) can be separated only when 
ul = cLr and ur — cLl- At the same time, by inspection 
of the complete system of four Dirac equations, 



(id 


- eA 


- 9^r 


- id r )u L y= 


muRy + 


d L Z, 

r 


(id 


-eA a 


+ gK 


+ id r )d[ J Z= 


mdRZ + 


u L — y, 

r 


(ido 


-eA 


-g^r 


+ id r )u R y= 


miiLy — 


dR Z, 

r 


(id 


-eA 


+ gK 


— id r )dftZ= 


mdhZ — 


ur — y, 

r 



where 

A± = (3 e ±-^U (C.8) 
one can see that the condition ul = dR and ur = d^ is 



inconsistent with the presence of the field H Q . The Dirac 
equation breaks up into two systems of equations for only 
two radial functions which are incompatible unless H r = 0. 

Nevertheless, just by inspection, one can see that the 
angular functions y kt m{8, and Z k ^ m (6,(p) that satisfy 
the equations, 



k-z k , m {e,ip) = -ky k , m (e,(p), 

^+y k , m {e, ^ = kz kim (e, v \ (c.9) 

do separate the angular variables in the Dirac equation 
(and do not separate them in Eq. (|C.5[0 . With this sep- 
aration of angular variables, Eqs. (|C.7|) yield a system of 
four differential equations for the four radial functions. 
Eqs. (|C.9I) are the equations for spherical harmonics but 
9 and ip are not the angles of the spatial angular coordi- 
nates. 
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